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Abstract 

In this paper, we consider the ID Navier-Stokes equations for viscous compressible 
and heat conducting fluids (i.e., the full Navier-Stokes equations). We get a unique 
global classical solution to the equations with large initial data and vacuum. Because of 
the strong nonlinearity and degeneration of the equations brought by the temperature 
equation and by vanishing of density (i.e., appearance of vacuum) respectively, to our 
best knowledge, there are only two results until now about global existence of solutions 
to the full Navier-Stokes equations with special pressure, viscosity and heat conductivity 
when vacuum appears (see [13| where the viscosity /i =const and the so-called variational 
solutions were obtained, and see [1] where the viscosity fi = fj,{p) degenerated when the 
density vanishes and the global weak solutions were got). It is open whether the global 
strong or classical solutions exist. By applying our ideas which were used in our former 
paper [S] to get 7?^— estimates of u and 9 (see Lemma [3.10[ Lemma [3.111 Lemma [3.121 
and the corresponding corollaries), we get the existence and uniqueness of the global 
classical solutions (see Theorem 11.11) . In fact, the existence of strong solutions would 
be done obviously by our estimates if the regularity of the initial data is assumed to be 
weaker. Like [8], we get iJ^— regularity of p and u (see Theorem ll.2p . We do not get 
further regularity of 9 such as i?''— regularity, because of the degeneration and strong 
nonlinearity brought by vacuum and the term {p,uux)x m the temperature equation. 
This can be viewed the first result on global classical solutions to the ID Navier-Stokes 
equations for viscous compressible and heat conducting fluids which may be large initial 
data and contain vacuum. 
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1 Introduction 

In this paper, we consider the Navier-Stokes equations for viscous compressible and heat 
conducting fluids (i.e. the full Navier-Stokes equations). The model, describing for instance 
the motion of gas, plays an important role in applied physics. Mathematically, the model in 
one dimension can be written as follows in sense of Eulerian coordinates: 

Pt + {pu)x = 0, p > 0, 

{pu)t + {pu^)x + Px = {pux)x, (1-1) 

{pE)t + {puE)x + {Pu)x = {pUUx)x + iK0x)x, 

for G (0, 1) X (0, +oo). Here p = p{x,t), u = u{x,t), P = P{p,9), E, 9 and k = k{p,9) 
denote the density, velocity, pressure, total energy, absolute temperature and coefficient of 
heat conduction, respectively. The total energy E = e + |n^, where e is the internal energy, 
/i > is the coefficient of viscosity. P and e satisfy the second principle of thermodynamics: 

In the present paper, we consider the initial and boundary conditions: 

{p, u, 9)\^^^ = {po, uo, 9o){x) in [0,1], (1.3) 

and 

{u, 9x)l^^, = 0, t>0. (1.4) 

Since the model is important, lots of works on the existence, uniqueness, regularity and 
asymptotic behavior of the solutions were done during the last five decades. While, because 
of the stronger nonlinearity in (jl.ip compared with the Navier-Stokes equations for isentropic 
fluids (no temperature equation), many known mathematical results mainly focused on the 
absence of vacuum (vacuum means p = 0), refer for instance to [T71 [TU [Ml ESI IMl EDI IM] for 
classical solutions. More precisely, the local classical solutions to the Navier-Stokes equations 
with heat-conducting fluid in Holder spaces was obtained respectively by Itaya in |17] for 
Cauchy problem and by Tani in [3l] for IBVP with inf po > 0, where the spatial dimension 
= 3. Using delicate energy methods in Sobolev spaces, Matsumura and Nishida showed in 
[291 [30] that the global classical solutions exist provided that the initial data is small in some 
sense and away from vacuum and the spatial dimension A'^ = 3. For large initial data and 
dimension = 1, Kazhikhov, Shelukhi in [25] (for polytropic perfect gas with p,K =const.) 
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and Kawohl in [24J (for real gas with k = k{p,9), /i = const.) respectively got global classical 
solutions to (ll.ip in Lagrangian coordinates with boundary condition ()1.4p and inf po > 0- 
The internal energy e and the coefficient of heat conduction k in [24J satisfy the following 
assumptions for p <'g and 6 >0 (we translate these conditions in Eulerian coordinates) 

eip, 0) > 0, 1/(1 + 6'-) < dee{p, 6) < iV(^)(l + 6^, 

< Ko{i + e'')<K{p,e)<Ki{i + e''), (1.5) 
^\dpK{p,e)\ + \dppKip,e)\ < ^(1 + ^"), 

where r £ [0,1], q > 2 + 2r, and v, N{'g), kq and ki are positive constants. For the perfect 
gas in the domain exterior to a ball in {N = 2 or 3) with fj,, k =const., Jiang in pTS] got 
the existence of global classical spherically symmetric large solutions in Holder spaces. 

In fact, Kawohl in |24] also considered the case of /i = ^(p) for another boundary condition 
with inf po > 0, where < p^ < p{p) < JIq for any p >0 and p^ and Jlq are positive constant. 
This result was generalized to the case p{p) = p" by Jiang in [19] for a G (0, j), and by Qin, 
Yao in [3T] for a G (0, ^), respectively. 

On the existence, asymptotic behavior of the weak solutions for full Navier-Stokes equa- 
tions (including the temperature equation) with inf po > 0, please refer for instance to 
[20\ [2T\ [23] for weak solutions in ID and for spherically symmetric weak solutions in bounded 
annular domains in {N = 2, 3), and refer to [12] for variational solutions in a bounded 
domain in (A^ = 2, 3). 

In the presence of vacuum (i.e. p may vanish), to our best knowledge, the mathematical 
results about global well-posedness of the full Navier-Stokes equations are usually limited to 
the existence of weak solutions with special pressure, viscosity and heat conductivity (see 
[U |13]). More precisely, Feireisl in [13] got the existence of so-called variational solutions 
in dimension N > 2. The temperature equation in [13] is satisfied only as an inequality. 
Anyway, this work in [13] is the very first attempt towards the existence of weak solutions to 
the full compressible Navier-Stokes equations in higher dimensions, where the viscosity p is 
constant and 
/■ 

K = k{0) E C2[0, oo), k{1 + 6") < K{e) < k{1 + 9") for all (9 > 0, 
p = p[p^ 0) = Ve{p) + 9Ve{p) for aU p > and 6* > 0, 
< Pe,7^e G C[0,oo)nCi(0,oo); Pe(0) = 0, Pe(0) = 0, (1-6) 
V'^ip) > aip~-^ - hi for ah p > 0; Ve{p) < a2p~ + h for all p > 0, 
Vg is non-decreasing in [0,oo); Ve{p) < 133(1 + p^) for all p > 0, 

where F < ^ if A'" = 2 and F = for A'' > 3; a > 2, 7 > 1, and ai, 02, 03, 61, k and 
K are positive constants. Note that the perfect gas equation of state (i.e. P = RpO for 
some constant i? > 0) is not involved in (jl.6p . In order that the equations are satisfied as 
equalities in the sense of distribution, Bresch and Desjardins in [I] proposed some different 
assumptions from [T3j , and obtained the existence of global weak solutions to the full Navier- 
Stokes equations with large initial data in or R^. In [l], the viscosity p = p(p) may vanish 
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when vacuum appears, and k, P and e are assumed to satisfy 

{K{p,e) = Ko{p,e){p + i){e^ + i), 

\p = rpe + p,{p), (1.7) 
l^e = C^O + edp), 

where a >2, r and C„ are two positive constants, Pc{p) = 0{p^^) and edp) = 0{p^^^^) (for 
some i > 1) when p is small enough, and Ko{p,6) is assumed to satisfy 

Co < Ko{p,0) < —, 
So 

for Cq > 0. We have to mention that the smooth solutions in ([0, oo); H'^{R^)) {d > 2+[^]) 
would blow up when the initial density is of nontrivial compact support (see |36j). On the 
local existence and uniqueness of strong solutions in M'^, please refer to for the perfect gas 
with ^, K =const. 

It is still open whether global strong (or classical) solutions exist when vacuum appears 
(i.e., the density may vanish). Our main concern here is to show the existence and uniqueness 
of global classical solutions to (|l.ip - p.4|) with vacuum and large initial data. In fact, the 
existence of the strong solutions to this problem is obvious if the regularity of initial data is 
assumed to be weaker. 

For compressible isentropic Navier-Stokes equations (i.e. no temperature equation), there 
are so many results about the well-posedness and asymptotic behaviors of the solutions 
when vacuum appears. Refer to [H [22l [261 [28] and [El [271 [351 [Ml [SHI [lQ| for global weak 
solutions with constant viscosity and with density-dependent viscosity, respectively. Refer 
to [6l [To] and [21 [3l [5l [33] for global strong solutions and for local strong (classical) solu- 
tions with constant viscosity, respectively. Recently, Huang, Li, Xin in [TU] and Ding, Wen, 
Yao, Zhu in [SI [9] independently got existence and uniqueness of global classical solutions, 
where the initial energy in [16] is assumed to be small in and p — p E C ([0, T]; i?^(M'^)) , 
u e C ([0,r];L»i(M3) nD3(]^3)^ P ^oo {\j-^T]-D'^{m?)) (for T > 0) which generalized the re- 
sults in [3], and the initial data in [8, 9J could be large for dimension = 1 and could be 
large but spherically symmetric for N >2, and {p,u) G C{[0,T]] H^{I)) {I is bounded in [8], 
and is bounded or an exterior domain in [9J). 

We would like to give some notations which will be used throughout the paper. 
Notations: 

(1) /= [0,1], a/ = {0,1}, Qt = Ix [0,T] forT>0. 

(2) For p G [l,oo], = LP{I) denotes the space with the norm || • \\lp. For k >1 and 
p G [l,oo], W'^^P = W'''P{I) denotes the Sobolev space, whose norm is denoted as || • 

ffk ^ p^fc,2(j)^ 

(3) For an integer A; > and < a < 1, let C^^'^{I) denote the Schauder space of 
functions on /, whose feth order derivative is Holder continuous with exponents a, with the 
norm || • ||c*;+c«. 
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In this paper, our assumptions are the following: 
(Ai): ljPo>0. 

{A2): IX = const. > 0, e = CoQ(0) + ec(/o), P = pQ{6) + Pc{p), k = k{6), for some constant 
Co>0. 

(A3): P^{p) > 0, edp) > 0, for p>0; Pee C2[0,oo); p\^\ < Ciedp), for some constant 
Ci > 0. 

(A4) : Q(-) e C2[0,oo) satisfies 

C2 (/3 + (1 - /3)e + 6^+') < Q{d) < C3 (/3 + (1 - P)6 + ^1+'^) , 
C4(l + 0") <Q'(e) <C5(1 + 0O, 
for some constants Cj > (i = 2, 3, 4, 5) and r > 0, /3 = or 1. 

(A5): K G C^[0,oo) satisfies 

C&{l + e'i)<K{e)<Cj{l + 6'i), 
for g > 2 + 2r, and some constants Cj > (z = 6, 7). 

{Aq): Q,Pc € C^[0,oo), and k satisfies 

\dlK{e)\<Cs{i + e'^-'), 

for > and some constant Cg > 0. 

Remark 1.1 (i) (Ai) is needed to get the upper bounds of 6 and 6t in terms of some norms 
by using mass conservation, Lemma \2.1\ and Corollary \2.1\ 

(ii) The case for the perfect gas (i.e. P = RpO, e = C^6 for constants R > and C^j > 0) 
is involved in the above assumptions. 

(Hi) As it mentioned in 124'}, the restriction on p (p=const., see other restrictions on k 
and e in i fi.5|) ) is not physically motivated. Physically, it seems more importantly that the 
state functions e, p and k usually depend on both p and 9. Particularly, the internal energy 
e grows as 9^^"^ with r ~ 0.5, the conductivity k grows as 9'^ with 4.5 < q < 5.5 and viscosity 
p increases like 9^ with 0.5 <p< 0.8 (see and references therein). Because of math- 

ematical technique, in the present paper, we assume p =const. and k = k{9) as in (see 
U.6\) ). From (A2)~(A^), we know that e and k grow respectively as 9^'^^ and 9'^, where q can 
be taken as q e [4.5, 5.5], and r can be taken as r = 0.5 if we consider 6 > 0. 

(iv) The restriction on q in (A^) (i.e. q > 2 + 2r) is same as il.5\) . and is the same as 
lll.6\) and (fi. 7p when we take r = 0. This assumption plays an important role in the analysis. 



Main results: 
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Theorem 1.1 In addition to {Ai)-{A5), we assume po > 0, po ^ , {^/f>o)x £ L°° , uq G 
n H^, 9q G H"^, dx9o\x=o,i = 0, and that the following compatible conditions are valid: 

pUOxx - [P{P0, Oo)]x = VPofl' /I o^ 

[i.o) 

[k,{0o)0Ox]x + P\u0x\'^ = 92, X £ I, 

for some 51,(72 £ L'^, cLnd (y^f^i)^ , (■y/po(72)^ G L^- Then for any T > there exists a 
unique global solution (p, u, 6) to hl.l\] - lT^\) such that 

p&C{%T]-H^), eC([0,r];/7i), ^gtyi,oo(g^)^ 

put GL~([0,r];ii'i), uteL\%T]-Hl), e L^{[0,T]- L^), 

peteL^{[0,T]-H^), 9eL^{[0,T]-H^), 9t e L^{[0,T]; H^). 

Remark 1.2 (i) il.8\) was proposed by Cho and Kim in to get local //^-regularity of u 
and 9 for the polytropic perfect gas. The detailed reasons why such conditions were needed 
can be found in Roughly speaking, gi and g2 are equivalent to yfput and i/pcf at t = 0, 
respectively. 

(a) By the Sobolev embedding theorems (cf. ^) and Lemma W^ we know from Theorem 

o 

p G C ([0, T]; n ([0, T]; , 

u G C([0,r];C2+-(/)) , {pu)t G C([0,r];C'^(/)), 

0GC([O,T];C2+-(/)), (pe)t GC([0,r];C"^(/)), 

for any T > and a G (0, |). This implies {p,u,9) is the classical solution to lll.l\) - [L4\ ). 

Theorem 1.2 In addition to {Ai)-{Aq), we assume po > 0, po ^ H^, {^/po)x G L°° , uq G 
H'^ n Hq, 9q G H^, dx9o\x=o,i = 0, q > 2 + 2r, and that the following compatible conditions 
are valid: 

PUoxx - [PiPO, Oo)]x = P093, 

(1.9) 

[l^{9o)9ox]^+ P\u0x\'^ = ^/P0 92, X G /, 

for some 53 G H^, (^/po^x-^s)^ ^ L^, and g2, [^92)^ e L^. Then for anyT > there exists 
a unique global solution (p, u, 9) to ^.1\) - [L4^ satisfying: 

pGC{[0,T];H^), pteC{[0,T];H'), ^gT^1'~(Qt), 
u G Ci[0,T];H^)nL^i[0,T]-H^), ut G i[0,T]; H^) D L\[0,T]; H^), 
{pu)t€C{[0,T];H^), ^pUxxt^L^{[^,T];L^), ^pet e L^{[Q,T]-L^), 
{pe)t G L-([0,T];i/i), 9 G L'^ {%T]- H^) rM^{%T]- H^), 9t G L\%T]-H^). 

Remark 1.3 (i) was proposed by Cho and Kim in where they consider the lo- 

cal existence of classical solutions for isentropic fluids (no temperature equation). Roughly 
speaking, g^ is equivalent to ut at t = 0. 



Global Classical Large Solutions to the Navier-Stokes Equations with Vacuum 



7 



(ii) We could not get 9 G C{[0,T]; H'^) (or L'^{[0,T]; H^)) even if [T^h is changed 
similarly to . because of the strong nonlinearity and degeneration brought by {puUx)x 

in the temperature equation and the appearance of vacuum, respectively. 

(Hi) Using ideas of Cho and Kim in 13], we can also get 

ueL^{[T,T];H'), eeL^{[T,T];H^), 

for T > 0. If we can obtain our estimates in higher dimensions, it will be useful to investigate 
the local (global) existence of classical solutions to the full Navier-Stokes equations (including 
the temperature equation) in M'^ (N > 2). For example, to guarantee jj.^p in (N = 2 
or 3) is valid for all t > 0, it is necessary to get 9 G L°° ([0, T]; H^^ . We will consider these 
problems in the near future. 

The constants Co in {A2) and the viscosity don't play any role in the analysis, we assume 
henceforth that Co = 1 and p = 1 for simplicity. 

The rest of the paper is organized as follows. In Section 2, we present some useful lemmas 
which will be used in the next sections. In Section 3, we prove Theorem [TTT] by giving the initial 
density and the initial temperature a lower bound 6 > 0, getting a sequence of approximate 
solutions to (ll.ip - (ll.4p . and taking 5 — t- 0^ after making some estimates uniformly for 5. 
More precisely, based on Lemma [2. II and the one-dimensional properties of the equations, we 
get /f^— estimates of the solutions. Using our ideas in [SI [9], we obtain //^—estimates of u 
and 9. In Section 4, using the similar arguments as in Section 3, we prove Theorem 11.21 

2 Preliminaries 

Lemma 2.1 Let Q = \a, b] be a bounded domain in M, and p be a non-negative function such 
that 

< M < / p<K, 
Jn 

for constants M > and K > 0. Then 



K 1 



M 



pv 



for any v G H^{Q). 



Proof. For any x G il, we have 



\v{x)\ < ^ 



'"(x) I p{y)dy 
n 



< 



M 



v{x)p{y)dy - / p{y)v{y)dy 



+ 



M 



p{y)v{y)dy 



< — 

- M 




nJy 



V({OdCp{y)dy 



1 



p{y)v{y)dy 



K 1 

< ^\\vx\\L^n) + 



M 



p{y)v{y)dy 



□ 
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Remark 2.1 The version of higher dimensions for Lemma \2.1\ can he found in fj^ or 113^ . 

Corollary 2.1 Consider the same conditions in Lemma \2.1\ and in addition assume 0, = I, 
and 

\\pv\\l^(i) < c. 

Then for any I > 0, there exists a positive constant C = C{M,K,l,c) such that 

\\v^\\l^{I) < CIK-y^ 11^2(7) +C, 

for any G H^{I). 

Proof. By Lemma l2.lt we have 

\\v^\\l^(i) < C\\{v^)^\\l2^i) + C j^p\v^\. 

Case 1: I G (0,1]. 

In this case, we use the Young inequality to get 



Il^h) < c\\{v%\\Lm) + c j^p\v\ + c j^p + c 

< C\\{v%\\l^j) + C. 



Case 2: / G (l,oo). 

In the case, we use the Young inequality again to get 



This gives 



HI) < C\\iv )x\\L^ii) + C\\v'-'\\L^^j) J^p\v\ 

< C\\{v%\\L2^j) + hv^\\L^^)+C. 



|f'||L-(/) < C\\{v^)^\\l2^i)+C. 



Lemma 2.2 For any v G Hq{I), we have 

lkl|L-(/) < I^xIIlI- 
Proof. Since w(0) = 0, we have for any x G I 



\v{x)\ = \v{x) — V 
This completes the proof. 



< ll^a;||Li(/)- 



Lemma 2.3 (\32^). Assume X C E C Y are Banach spaces and X 
following imbedding are compact: 



□ 



□ 



E. Then the 



(i) |99:<^GL''(0,r;X),^ GLi(0,T;y)| -^-^L'?(0,r;i?), if 1 < g < oo; 



{ii) : V e L'^{0,T;X),-^ e L'{0,T;Y)} C{[0,T];E), if 1< r < oo. 



Global Classical Large Solutions to the Navier-Stokes Equations with Vacuum 



9 



3 



Proof of Theorem 11.1 



In this section, we get a global solution to (jl.lj) - (jl.4l) with initial density and initial tem- 
perature having a respectively lower bound (5 > by using some a priori estimates of the 
solutions based on the local existence. Theorem 11.11 would be got after we make some a priori 
estimates uniformly for 6 and take (5 — )• 0^. 

Denote pfj = po + 5 and 9^ = 6o + 6 fov 6 £ (0, 1). Throughout this section, we denote 
c to be a generic constant depending on po, uq, ^O; T and some other known constants but 
independent of S for any 5 G (0, 1). 

Before proving Theorem ll.il we need the following auxiliary theorem. 

Theorem 3.1 Consider the same assumptions as in Theorem \l.l[ Then for any T > and 
S € (0, 1) there exists a unique global solution {p,u,9) to il.l\ )- [T^\ ) with initial data replaced 
by (pq,uo,0()), such that 



peC{[0,T];H^), pteC{[0,T];H^), pu e L^^T]; L^), p > > 0, 
u E C{[Q,T];H^r\Hl), ut e C{[^,T];H^)r\L\[d,T]-H^), uu G L2([0, T]; L^), 
^>Q>0, 0GC([O,r];i/3), gC'([0,T];Fi)nL2([0,T];ii'2), Ou & L^{[Q,T\-L^), 



where cs is a constant depending on 6, but independent of u. 

Proof of Theorem I3.lt 

The local solutions as in Theorem 13.11 can be obtained by the successive approximations 
like in [3]. We omit it here for brevity. The regularities guarantee the uniqueness (refer for 
instance to [3]). Based on it. Theorem 13 . 1 1 can be proved by some a priori estimates globally 
in time. 

For any given T G (0,+oo), let {p,u,6) be the solution to (|l.ip - (|1.4p as in Theorem 13.11 
Then we have the following basic energy estimate. 

Lemma 3.1 Under the conditions of Theorem \3.1{ it holds for any < t <T 



Proof. Integrating (jl.ip i and (jl.iP '^ over /x [0, t], and using ()1.4p , {A2) and {A4), we complete 



Proof. The proof of the upper bound of p relies on constant viscosity (i.e. p = const.). It is 
similar to [37] . 




the proof of Lemma |3. 11 



□ 



Lemma 3.2 Under the conditions of Theorem \3.1l it holds for any {x,t) £ Qt 



< p{x,t) < c, 
e{x,t) > 0. 



Denote 




(3.1) 
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Differentiating (j3.ip with respect to x, and using (jl.ip -;. we have 

Wx = pu. 

This together with Lemma 13.11 and the Cauchy inequality gives 

\Wx\ < C. 

It follows from ([H]), ([O]), (^2), (^3), (^4), dLlI), and Lemma O that 

w < c. 

J I 

This gives for any (x, t) G Qt 



w{x,t)\ < 


w{x, t) — J w 


+ 













< 



+ c 



< J \Wx\ + C < C, 

which implies 

For any {x,t) £ Qt, let X{s;x,t) satisfy 

^^^^i^ = u{X{s;x,t),s), 0<s<t, 
X(t; x, t) = X. 

Denote 

F{x, t) = exp {w{x, t)} 

It is easy to verify 



djpF) iX{s;x,t),s) 
ds 

Multiplying (j3.4p by exp (Jq P), we have 

d 



/ dp dw 
-pPF. 



0. 



Integrating it over (0,t), we have 
which implies 

p{x,t) > 0, 

for any {x,t) £ Qt- 

By (j3.2p . (j3.5p and P > 0, we get the upper bound of p. The lower bound of 1 
by (j3.7p and the maximum principle for parabolic equations. 
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Lemma 3.3 Under the conditions of Theorem \3.1[ it holds for any given a G (0, 1) 



'T 



where c may depend on a. 



Remark 3.1 a was usually taken as 1 when the basic energy inequality was done (see fl]/ 
and references therein). This depends on pologOo £ which can not be got under the 
assumptions of Theorem \1.1\ and Theorem since 6q may vanish. 

Proof. From ()1.2p and (jl.ip . we get 

peeOt + pueee^ + OPeu^ = ul + {^{6)9^)^ . (3.6) 
Substituting e = Q{e) + edp) and P = pQ{e) + Pc{p) into ([32]), we get 

pQ'{e)et + puQ'{e)e, + peQ'{e)u, = ul + {K{e)e^)^ , (3.7) 



or 



{pQ)t + (puQ),, + peQ'{e)u, = 4 + {K{e)e,)^ . (3.8) 



Multiplying ()3.7p hy 9 integrating the resulting equation over Qt, and using integration 
by parts, we have 



P 



I 



Q'iO 



cc 
'5 



Ji Jo ? JQt 



T 



< c I p{i + e'^n + c + l [ $ + cl 



/ ^ JQt P JQ. 



T 



1 / UZ 



2 rT 



< c+- -^+c maxe^+'-", (3.9) 

where we have used (A4), the Cauchy inequality, Lemma 13.11 and Lemma 13.21 Now we 
estimate the last term of (13.91) as follows: 



c / max0i+'-" < c+ / \\9''~''94l2 
'0 ^e-f Jo 



< c + c 



T / f f)2g2r-a+l \ \ 







2 



2]q^ 91+" ^ ' 

where we have used Corollarv 12.11 Lemma |3.H (A5) and the Cauchy inequality. By (|3.9p . 
p.lOP and (^45), we complete the proof. □ 



Global Classical Large Solutions to the Navier-Stokes Equations with Vacuum 12 
Corollary 3.1 Under the conditions of Theorem \3.1\ it holds 





Proof. By Corollary 12.11 and Lemma |3. 11 we have 

i-T t-T 



j 

Jo 







< c 



T 



Ji 

T 



q — a — \ 

" Ox 



2 



J I 



c / / o'^-'^'^el + c 



< c. 



□ 



Lemma 3.4 Under the conditions of Theorem Vj.l\ it holds 



ul < c. 

Qt 

Proof. From (jl.ip i and ()l.ip 9. we get 

put + pUUx + Px = Uxx- (3-11) 

Multiplying (|3.11|) by u, integrating it over /, and using integration by parts, we have 



< ^ i< + c J^p^Q' + c J^P^ 



< - jul + c J^e^+^^ + c, 



1 

2 

where we have used the Cauchy inequality, (^2)5 (^s); (^4) and Lemma |3.2[ This implies 

d 



, pu^+ < <csupe^-"+^ + c. 
J I J I xei 

Integrating it over (0,t), and using Corollarv 13.11 we complete the proof of Lemma 3.4. □ 
Lemma 3.5 Under the conditions of Theorem \3 . li it holds for any <t <T 

[ {ul + pe'^^'^n + [ {pul + {i + e'^fel)<c. 

J I JQt 

Proof. Multiplying p. lip by ut, integrating it over /, and using integration by parts, Lemma 
7I[ Lemma 13.21 and the Cauchy inequality, we have 



= ij/^^^- j^Pnn^^t- jPtUx 



1/" 2,1/" 22 



< 
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which implies 

We are going to estimate the last term of the right side of (|3.12|) . Using (^2), (I3.8|) . (Il.ip i 
and integration by parts, we have 

-2j^Ptiu,-P) = -2 j{pQ)t{u^-P)-2 j{P,)t{u^-P) 

= -2 [{kO^)^ + ul- (puQ)^ - p9Q'{e)u^] (n^ - P) 

+2 j ^P'^{p){p^u + pu^){u^ - P) 

= j l^^xiUxx - Px) -2 j ul{ux - P) -2 J puQ{uxx - Px) 

+2 ^ peQ'{e)Ux{Ux - P) - 2 ^ PcU{Uxx -Px)-2j^ PcUxiUx - P) 

+2 j^pP',{p)ux{ux- P). 

This, together with ()3.1ip . (^2), (^4); Lemma [221 Lemma [321 the Cauchy inequality, and 
W^'^{I) ^ L'^il), gives 

-'^j Pt{ux-P) - i^(^x{put + puux) + 2\\ux - P\\l°° j ul,-2 j puQ{put + puux) 

+csup(l + 0^+'^) / ul + csnv{l + d^+'') I pQ"^ + csupe^+'' 
x&i J I xei J I xei 

—2 J Pcu{put + puUx) + c J ul + c J pQ^ + c 

- \j^P^^t+cj^K^Ol + c(^jul^ +c{\\Ux- P\\L^ + \\pUt + pUUxWl^) j^ul 

+c [ ul [ pQ'^ + csup{l + 9^+"^) [{ul + pQ^) + csup9^+'' + c 
Ji Ji xei Ji x&i 



jul^ +^ j^pul + c jul j^pQ^ 
+csup(l + 0^+'^) [iul + pQ^) + csupe^-^'' + c. (3.13) 

x£l J I XGl 



Substituting (|3.13|) into (|3.12p . we have 



+csup(l + ^1+'') / {ul + + csup^i+'' + c. (3.14) 
xei Ji xei 

Integrating (|3.14p over (0,t), and using {A2)-{A4), Lemma [3^ Corollarv l3.1l Lemma [3^ and 



Global Classical Large Solutions to the Navier-Stokes Equations with Vacuum 



14 



the Cauchy inequality, we have 



1 
2 

< c 



Ji 



PUt + / 
2 



fiful) +2 [{pQ + P^)n^ + c f I n^el + c f ! ul I pO^^^^ 
Jo \Ji / JI Jo Ji Jo Jl JI 



+c / sup 6'^+'' I ul + c I sup(l + 6' 



xel 



^l+r■^ 



xGl 



+ C 



< C 



32+2r 



+ C 




Jl 



K^el + c 




Jl 



ul / p9 



32+2r 



Jo x£i Jl Jo xei Jl 



+ c. 



The second term of the right side can be absorbed by the left. After that, we have 

r-t 




Jl 



< c 



K^el + c 



Ji 




sup 01+'^ f ul + c f SUp(l+^l+'^) / p^2+2r^^^ 
x&I Jl Jo x£l Jl 



32+2r 



(3.15) 



Here, we have used Lemma 13.21 and the Young inequality on the second term of the right 
side. Note that the terms about 6 in (|3.15p need to be handled. To do this, we make use of 
(I3TD. 

Multiplying (|3.7p by K,iOdS,, integrating it over /, and using integration by parts, (^44) 
and (^5), we have 



dtl " 



Q'iv) / ^iOd^drj 



+ / K 



,2/)2 



ul / KiOd^ - / p9Q'{9)u^ / Ki^dC 



<c\\(l + 



ei)9U^ ^ ul + c||(l + 9'^)e\\Loo J^p{l + 9'+')\u,\. 



By Corollary O and (^5), we get 



+ 9'^)9\\loo < c\\Ke,\\L2 + c. 



(3.16) 



(3.17) 



Substituting ()3.17p into (|3.16p . and using the Holder inequality, the Cauchy inequality and 
Lemma 13.21 we get 



d_ 

di 



+ / K 



Q'iv) / md^dr] 
Jo 

< c\\k9,Al^ jul + cjul + c\\k9,\\l2 J^p{l + 9^+'-)\u,\ + c^p(l + 9^+^)\u^\ 



< c\\k9,\\l2( I ul + \\pil + 9'+'-)h2\\u,h2]+c I ut + c I pil + 9'+'n 

2 



1 



< =; I K^Ol + c ( / + c 



32+2r 



ul + c 



-,2+2r 
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which implies 



d_ 

dt 



P 



Jo 



2 

< c' ' 



Integrating it over (0,t), and using {A4), (^5), Lemma l3.ll and Corollary 13. H we get 

By (|3.15|) , (|3.18|) , Corollary 13. H Lemma \'6A\ and the Gronwall inequality, we complete the 
proof. □ 

Lemma 3.6 Under the conditions of Theorem \3.1{ it holds for any < t <T 

[ {pi + Pt) + [ ulx< c. 
Jl JQt 

Proof. Differentiating (jl.ip i with respect to x, we have 

Pxt + PxxU + 2pxUx+ pUa;a: = Q. (3.19) 

Multiplying (|3.19p by 2px^ integrating it over / and using integration by parts, we have 



d 

'dt J I 



-T. I pI = -3 / pIux - 2 / ppxu^ 



= -3 J pI{ux - P) -3 j pIp -2 j ppxUxx 

< c{\\ux-P\\l2 + \\uxx - PxWl^) j^pl + c j^ulx + c j^pl 

< c{l + \\^put\\L2) j^pl + c julx, (3.20) 

where we have used (|3.1ip . the Sobolev inequality, (^2)-(^4)) the Cauchy inequality. Lemma 
12.21 Lemma 13.21 and Lemma 13.51 

It follows from (j3.1ip . Lemma 12.21 Lemma 13.21 [A2)-{A/i)., Lemma 13.51 and the Cauchy 
inequality that 

julx < cj^pu't+cl^l^ul^ +c j^plQ^ + c j^p'[Q'{e)fel + c j^pl + c 

< C [ pul + CSUp(l + 02+2r) /■ ^2 ^ , /■(! + e^lfol + c 

Jl xel Jl Jl 

< c [ puj + csup(l + 09-"+i) [ pI + c [{1 + e'i fel + c. (3.21) 

J I xei Ji Ji 

Substituting ()3.2ip into ()3.20p . and using the Gronwall inequality, Corollarv 13.11 and Lemma 
13.51 we get 

%'<c. (3.22) 
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By IK2n . IK22h . Corollary O and Lemma [331 we have 

It follows from (jl.ip i . (j3.22p . Lemma 12.21 Lemma 13.21 and Lemma 13.51 that 

Pt < c. 

The proof of the lemma is complete. □ 
Lemma 3.7 Under the conditions of Theorem \3.1l it holds for any < t <T 



J I JQt 
Proof. Differentiating (13. lip with respect to t, we have 

PUtt + PtUt + PtUUx + PUtUx + pUUxt + Pxt = Uxxt- (3.23) 

Multiplying (|3.23p by ut, integrating the resulting equation over /, we have 



2di 



oil I PUt+ ult 



-2 / pUUtUxt - / PtUUxUt - I pufux + / PtUtx 



PUt 



< 2\\^Ut\\L2\\^u\\L^\\Uxt\\L2 + \\ut\\L°°\\u\\L-^\\Pt\\L2\\Ux\\L2 + \\u^^^ j I 
+ \\Pcip)\\L°°\\Pt\\L2\\Uxt\\L'2 + \\Q{d)\\L^\\Pt\\L2\\Uxt\\L'2 + l|pQ'(6')6't ||l2 ll^xt |Il2 

< \(ul, + c( puUc + c( [ pu\ + csup02+2r + e f p (l + 0^+'-) Q\. 

^ J I J I Jl J I xel J I 

Here, we have used p.ip i . integration by parts, the Holder inequality, the Cauchy inequality, 
the Sobolev inequality, {A2)-{Ai), Lemma 12.21 Lemma 13.21 Lemma 13.51 and Lemma 13.61 
The first term of the right side can be absorbed by the left. This implies 

< c f puj + c + c f f puf + csup ^2+2^ + c / p (1 + 0"+'^) el (3.24) 

Jl Jl J I xel Jl 

Integrating (|3.24p over (0, t), and using Corollarv 13.11 and Lemma 13.51 we have 

/ pn? + f [ul< [ pu'tiO) + C + C f [uiJ pu^t+c f [ p{l + 05+0 (3.25) 
Ji Jo Ji Ji Jo Jl Jl Jo Jl 

Multiplying ()3.1ip by taking t — )- 0+ and using (jl.SP i . we have 

I^Oxx — 'P{Poj^o)x\ I : 

\^Ut{x,{))\ < J = + \J Po\uoUOx\ 



, \uoxx - P{po,%)x\ , \P{po,(^o)x - P{Po,Go)x\ , r~s\ I 
< — \ 1= — ^^—^ \- y Po\uouox\ 



< \gi\ + c^={\pox\ + \dox\) + c, 
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which implies 

pufiO) < c. (3.26) 



Substituting (3.26) into ()3.25p . we have 

/ pn? + f [ul<c + c f [ uL / pu? + c /* / p (1 + 0^+0 el (3.27) 
Ji Jo Ji Jo Ji Ji Jo Ji 



Multiplying (|3.7|) by |^Jq i^iO^^j (i-^- '^(^)^O) integrating the resulting equation over /, and 
using integration by parts, (^44), (^45), Lemma 12.21 Lemma 13.21 Lemma 13.51 and the Cauchy 
inequality, we have for any e > 



< ^J^ pQ'K0f + cj^ pu^Q'kOI + c pQ'^Q'kuI 

< II pQ'^Ot +cj{i + e''f6l + c(i + j^ uL) p{i + 0''"'^') 

Ji Jo J Ji xei 



d_ 



which combining Lemma 12.21 Lemma 13.21 Lemma 13.51 implies 

j^pQ>{e)K{6)d^, + jJ^KHl 



< C 

Integrating it over (0, t), and using (^4) and (^5), Lemma 12.21 Lemma 13.51 Lemma 13.61 and 
the Cauchy inequality, we obtain 



< c 



Jo Ji Ji 

I u1 / K(^)d^ + ce / / uli + Cs 
I Jo Jo Ji 



rt 

x£l Jo Jl 



< csup{l + e'^)e + ce I I ul^ + Ce 



t 

2 



< c\\{i + e'i)e^\\L2 + ce jui^ + Ce 



After the first term of the right side is absorbed by the left, we get 



< 
- 2 



f [p{l + + /(I + <ce f [ ul, + c (3.28) 

JO Ji Ji Jo Ji 
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Multiplying (j3.28p by 2c, adding the resulting inequality to (j3.25p . taking e = and using 
the Gronwall inequality and Lemma 13. 6^ we complete the proof of Lemma 13.71 □ 

Prom Corollary 12.11 Lemma 13.11 and Lemma 13. 7^ we get the following corollary immedi- 
ately. 

Corollary 3.2 Under the conditions of Theorem \3.1\ it holds 

P\\l°-{Qt) < c. 

Corollary 3.3 Under the conditions of Theorem \3.1\ it holds for any <t <T 

\\u\\wi,^(Qt) + / ul^ + el^< c. 

J I JQt 

Proof. It follows from (|3.2ip . Lemma 13.61 Lemma 13.71 and Corollary 13.21 that 

which, combining Lemma l2.2( Lemma 13.51 and the Sobolev inequality, gives 

\W\\W'^{Qt) ^ c. (3.29) 

By (I32D, Corollary [321 (^4), (^5), Lemma US Lemma [321 dSSHj), Lemma [321 the Holder 
inequality, Sobolev inequality and Cauchy inequality, we have 



/C < c f ei + c f ui+ f pel + c f uHl + c f i 
ji Ji Ji Ji Ji ji 

< Cp^O^^L^ I OI+ [ P0t+C 



I J I 

< c\\9^^\\l2 + / pOf + c 



After the first term of the right side is absorbed by the left, we get 



4 < / pd^ + c. (3.30) 
/ J I 



Integrating (j3.30p over [0,T], and using Lemma 13.71 we get 



oL < c. 

'Qt 

This proves Corollary 3.3. □ 
Lemma 3.8 Under the conditions of Theorem \3.1{ it holds for any <t <T 

\\p\\wi'°°iQT) + \\Pt\\L°°(QT) + / (pIx + pit) + / (Ptt + '^Ixx) ^ C. 

JI JQt 
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Proof. Differentiating (j3.19p w.r.t. x, we have 

Pxxt — PxxxU 'ipxx^x '^Px'^xx PUxxx- (3.31) 

Multiplying (j3.3ip by 2pxx, integrating it over /, and using integration by parts and the 
Holder inequality, we have 



dt 



I 



-Jl I Pix = -5 / pi^Ux - 6 / PxPxxUxx - 2 / pPxxUx 



< b\\Ux\\L°° I pIx + ^\\Px\\l°°\\Pxx\\l2\\Uxx\\l'^ + '2\\p\\l°°\\Pxx\\l'^\\Uxxx\\l'^- 



By the Sobolev inequality, Cauchy inequality. Lemma [321 Lemma ESI and Corollary 13.31 we 
have 

j^plx<c j^pl^ + c jul^^ + c. (3.32) 

The next step is to estimate the term JiU^^x- Differentiating (|3.11|) with respect to x, we 
have 

Uxxx = PxUt + pUxt + PxUUx + pul + pUUxx + {Pc)xx + {pQ)xx- (3.33) 

By (^3), (^4), Lemma [2.21 Lemma [3.21 Lemma [3.61 Corollarv 13.21 Corollarv 13.31 and the 
Sobolev inequality, we get 

ulxx < c / p^ul^ + c / pluj + c pI^ + c el^ + c 



I JI JI JI JI 

2 , „ / „2 , „ a2 



< c jut^ + c pi^ + c + c. (3.34) 

Substituting (|3.34p into (|3.32p , and using the Gronwall inequality. Lemma 13.71 and Corollary 
we get 

%L<c. (3.35) 



By (|3.35p , Lemma 13.21 Lemma 13.61 and the Sobolev inequality, we have 

\\p\\w^'^{Qt) ^ c. (3.36) 
By ([CTp . (1335|), Lemma [321 and Corollary [331 we get 

/ ulxx < C. 

JQt 

The estimates of pxt and pu can be obtained directly by (j3.19p . (jl.ip i . (|3.35p . (|3.36p . Lemma 
12.21 Lemma 13.21 Lemma 13.61 Lemma 13.71 and Corollary 13.31 The proof of Lemma 13.81 is 
complete. □ 

Lemma 3.9 Under the conditions of Theorem \3.1{ it holds for any <t <T 

[ pel+ [ \{Kex)tf<c. 

J I JQt 
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Differentiating (j3.7p w.r.t. t, we have 

pQ'ett + pQ"el + ptQ'dt + {puQ'9^)t + {peQ'u^)t = lu^u^t + {^.^xU- (3.37) 

Multiplying (lOTl) by (/g K{S,)d(,)t (i.e. ^{6)61), integrating it over /, and using integration by 
parts, p.ip i ■ (^4), (^5), Corollary 13.21 Lemma [3?2t CoroUarv 13.31 and the Holder inequality, 
we have 



ld_ 

2di 



|2 
)t\ 



= ^ ptQ'ne^t pQ"efK + i ^ pQ'K'ef - j{puQ'e^)tKet 

- j {p9Q'ux)tK9t + 2 j UxUxt^Ot 
< \ J^{pu),Q'k9^ + cWKOthoo j^pO^t - J^puQ'{K0M - J^puiQ"K - Q'k')9^0, 

\ {pu)tQ'9xH9t + C ult+C p9f - / Pt9Q'UxK9t + c\\K9t\\L^\\U:„t\\L2\\Ux\ 

I Jl Jl Jl 



This, combining integration by parts, (^4), (^5), Lemma 12. H Lemma 12. 2| Corollary 
Corollarv 13.31 Lemma 13.71 Lemma 13.81 and the Cauchy inequality, gives 



2di 



)t\ 



c\\m\\L^ I^p9t 



j^pQ'K9l + j^\{K9^ 

< - ^ ^ puQ"e^K9l -^j^ puQ'k'9^9^ - ^ puQ'K9t9^t + 

+c\\^/p9t\\L4{'^0x)t\\L^ + c\\9Jl^ j^p9^ + c\\K9t\\L^ + C j^ult + C j^p9^ 
+c\\K9t\\L°°\\Uxt\\L'2 

< c\\9x4l^ j^p9f + c\\K9t\\L^ j^p9j + c\\^9t\\LA\{^^Gx)t\\L^ + c\\K9t\\L^ 

+C jult+C j^p9j + c\\K9t\\L'^\\Uxt\\L2 

< c\%x\\l^ j^p9''^ + c {\\{K9t)x\\L^ + I^P^\Ot\^ j^pOl + c\\^p9t\\LA\{^^0,)t\\L2 



+c||(At^t)x||L2 +C j^pK\9t\+C juit + C j^p9t+c[\\{K9t)x\\L^+ J^P'^\^t\j hxth^, 

which together with the Cauchy inequality. Lemma 13.21 (^5) and Corollarv 13.21 gives 



where we have used {K9t)x = {i^9x)t- This gives 



2 

2 I „ / / „fl2 \ , „ / „.2 



+ c p9t\ +c u%t + c, 



pQ'k91 + I^\{k9,),\^ <cl^9l, + c^p9^^ +cl^ul, + c. 



Integrating it over (0,t), and using (A4), (^5), Lemma 13.71 and Corollarv 13.31 we obtain 

_^/5^? + ^*y^l(^^x)d'<c^P^'(0) + cyJ|*(^y^pe?) +c. (3.38) 
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Multiplying (13.70 by , taking t — > 0+, and using (|1.8p 9. we have 



which implies 



< c\g2\ + ^={l + \eo^x\) + c, 

^0 



< c. (3.39) 

Substituting (j3.39|) into (|3.38|) . using the Gronwall inequality and Lemma [221 we complete 
the proof. □ 

Corollary 3.4 Under the conditions of Theorem \3.1\ it holds 

rT 

IMlo. < C. 







Proof. By Lemma l3.2t (A^), Corollary 12.11 Corollary Lemma |3.9|, and {K9t)x = {K9x)t, 
we get 

/ IlK^tllioo < c / \\{Ket)x\\l2 + c < c. 

Jo Jo 

This combining (^5) completes the proof. □ 
Corollary 3.5 Under the conditions of Theorem \3.1\ it holds 



Proof. Since 



we obtain 



KOxt — {l^dx)t — KOtOx, 



[ < c f 



< c / {{kOxU^ + c / {k') 



t "x 

T 

2 / a1 

x 



< C + C sup 

Jo xei Ji 

< c, 



where we have used (A5), Lemma [5771 Lemma [331 Corollary 13.21 and Corollary 13.41 □ 
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Corollary 3.6 Under the conditions of Theorem \3.1\ it holds for any <t <T 

\\^\\w'^-°°(Qt) + / + / ^Lx < C. 
Jl JQt 

Proof. From (|3.30p and Lemma 13.91 we have 

< c, (3.40) 
which, combining Corollary 13. 2t Lemma 13.71 and the Sobolev inequality, gives 

I|6'||wi.-(Qt) ^ c. (3.41) 

Differentiating (|3.7p w.r.t. x, we have 



i^'dxxx — —'^k'OxOxx — k"6^ — 2uxUxx + pQ'Gxt + PxQ' Gt + pQ"Gx(^t 

+{puQ'0x)x + {pOQ'ux)x. (3.42) 

By dS^Ql), dSai]), ([332]), (A4), (A), Lemma [33 Lemma [3J] and Corollary [331 we have 

ixx < c [ elel, + c [el + c [ ulul, + c / p^el + c / /^^^^ + , f ^2^2^2 
7 J/ J/ Ji Ji Ji Ji 



+c j \ ipuQ'9x)xf + c j \ {p9Q'ux)x^ + c 



< c I pHlt + c / ple^t + c 



< c el^ + csupe^ + c. (3.43) 
By ()3.43p . CoroUarv 13.41 and Corollarv 13.51 we obtain 

□ 



The next lemma, which we used in [8] to get iif^— estimates of velocity, plays an important 
role in getting estimates of 9 in the following. 

Lemma 3.10 Under the conditions of Theorem \3.1l it holds 

\\(.Vp)x\\l^(Qt) + \\{y/p)t\\L°-{QT) < C. 
Proof. Multiplying (ILlh i by , we have 

^Vp 

{Vp)t + iVp)xU + ^^/pUx = 0. (3.44) 

Differentiating ()3.44p with respect to x, we get 

3 1 

{\^)xt + {yfp)xxU + ■^{y/p)xUx + -^y/pUxx = 0. 
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Denote h = {y/p)x, we have 

3 1 

ht + h^u + -hUx + -^yfpUxx = 0, 

which implies 

^|/iexpQ^ dxu{X{T;x,t),T) dr^^ = -^^{dxu)exp j dxu{X{T;x,t),T) drj , 

(3.45) 

where X{s]x,t) is the solution to (|3.3p . 
Integrating (|3.45|) over (0,t), we get 

h{x,t) = exp ^-^ J dxu{X{T;x,t),T) dr^ h {X{0; x,t),0) 

-^exp^-^y dxu{X{T;x,t),T)dT^ j ^{dxu)exp(^j dxu{X{T;x,t),T) dr 

This together with Corollary 13.31 Lemma 13.81 and the Sobolev inequality, implies 

UVPUL^{Qr)<(^- (3-46) 
From ([3^ . Lemma [32] and Corollary [331 we get 

\\{Vp)t\\L°-(QT) < C. 

This proves Lemma 3.10. □ 

The next lemma will be used to get if^— estimates of 9. 
Lemma 3.11 Under the conditions of Theorem \3.1l it holds for any <t <T 



I 



Proof. Multiply (I3.37P by p^'^{K,6t)t (i-e. p'~'^ kOu + p^^k'O^, where 71 is to be decided later), 
and using integration by parts, we have 



2 



+2 j uxuxtip^' kOu + P^^K'e^t) - J^p"''+'Q'K'6^ett 

- y {pQ"Ot + PtQ'Ot + {puQ'9x)t + {p9Q'ux)t) {p^' kOu + p"" k9^) . (3.47) 

We are going to look for the minimal of 71. It seems that the second term of the right side 
plays an important role. 



71 / P^^ ^PxK'9u{n9x)t = -271 j p^^ ^{y/p)xt<'9u{K,9x)t 

< \j^p^-+^KQ'9l + cj^p^--\Kex)t\\ (3.48) 
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where we have used Lemma [3. 101 (^4), (^5), Corollary 13.21 and the Cauchy inequality. 

From Lemma 13.91 we know that Jq^ \ {n9x)t\'^ < c. This implies that the minimal of 71 
should be 2. Substituting 71 = 2 into (|3.47p and (|3.48p . and then substituting ()3.48p into 
([OTl) . we have 



1 ( 1 

2 / /■ A 2 



+ci^j^p^KQ'9lj |1+ iyJ^P^tj +\\Oxt\\L^ + \\VpUt\\L^ + \M\L^ + \\VpBt\\L^ ^ 

where we have used (A4), (A5), Lemma [3771 Lemma [3^ Lemma [3^ Corollarv 13.31 Corollary 
13.61 and the Cauchy inequality. This implies 

Integrating it over (0, t), and using (^4), (^5), Lemma ISTTl Lemma [3^ Corollary 13.41 and 
Corollarv 13.51 we have 



p'\{^Ox)t\" + / * / P^el < c I p'\{K6x)t\' (0) + c. (3.49) 
Jo Ji Ji 

By (yls), ([339]) . ([3:i2]l and Corollary [321 we get 

p'\{nex)t\'{Q) < cJ^p'elM + cJ^pe^io) 

< ^1100111-3 + 01114011^2+0 J pluf{0) + c 

< c + cJ^li^U^puUo) 

< c. (3.50) 
Substituting (j3.50p into (j3.49p . we complete the proof. □ 
Corollary 3.7 Under the conditions of Theorem \3.1\ it holds for any <t <T 

l{oL. + p'olt)<c. 

Proof. A direct calculation gives 

pKO^t = p{i^Ox)t - pn'OtOx, 

which implies 

I^P^ol < cJ^p^\i^e,)t\' + c\\e,\\lo. I^pe^ 

< c. (3.51) 
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Here we have used (^5), Lemma 13.81 Lemma 13.91 Lemma 13.111 and Corollary [37 
From the second inequality of (|3.43p . we obtain 

< c, 

where we have used Lemma |3.9[ (|3.5ip and Lemma |3.10[ □ 

The next lemma will be used to get i/^— estimates of u. 
Lemma 3.12 Under the conditions of Theorem \3.1[ it holds for any <t <T 

JQt 

Proof. Similarly to Lemma I3.1H multiplying (j3.23p by p^uu, and integrating it over J, we 
have 

= j pptulf- -2 j pPxU^tUtt - j P^UttiptUt + PtUUx + pUtUx + pUUxt + Pxt) 

< C jult-A j^p^ {y/p)xUxtUtt + P^ul +cj^pU^+cj^\ ipQ)xt\'^ + C j ^\{Pc) xt\^ + C 

- \j]p^''tt + cjult + c\\et\\lo. + cj^elt + c. 

Here, we have used integration by parts, Lemma [3. 7| Lemma [3.81 Lemma [3.10l Corollary 13.31 
Corollary 13.61 and the Cauchy inequality. 

The first term of the right side can be absorbed by the left. After that, we have 

.,3„,2 I ^ /" ^2,2 ^ ^ f „,2 I i\n ||2 I „ /" o2 



^ P + ^ J^P u^t<c + c||6'i ll^oo +cj^ + c. 

Integrating this inequality on both side over (0,t), and using Lemma [221 Corollary 13.41 and 
Corollary 13. 5[ we have 

/* / p\l + [ p^l, < [ p^lM + c. (3.52) 
Jo Ji Ji Ji 

Similarly to ([330)1 . we use ([3:26]) . ([3331) . (^3) and (^14) to get 

j^P^ultiO) < c\\uofH3+c\\eo\\l2+c\\pofH2+cJ^puf{0)+c 

< c. (3.53) 

Substituting (j3.53p into (j3.52p . we complete the proof. □ 

By (j3.34p , Lemma 13.71 Lemma 13.81 Lemma 13.101 Lemma 13.121 and Corollary 13.61 we get 
the following corollary. 
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Corollary 3.8 Under the conditions of Theorem \3.1\ it holds for any <t <T 



11^ < r 

"'XXX — ^' 



Prom the above estimates, we get 

\\{Vp)x\\l'^ + + WpWh^ + 11/0*11//! + \Mh3 + llmllz/i + WVP'^th^ + \Mh3 

+ \\VP0t\\L^ + WpGtWm + [ Ht + pI + 0^ + Olt + P'^tt + P'f^tt) < c. (3.54) 

JQt 

Corollary 3.9 Under the conditions of Theorem \3.1{ there exists a positive constant cs de- 
pending on 6 such that for any {x,t) G Qt, it holds 

p(x,t) > - > 0, 

' ^ - c (3.55) 

0{x,t) >cs>0. 

Proof. By (13. 5|) . (^3), (A^), Lemma [3^ and Corollary 13. 6^ we have for any {x,t) G Qt 

p{x,t) > -. 

c 

This gets p.55P i . (I3.55P 9 can be got by (|3.55P i . (j3.54p . ()3.7p and the maximum principle for 
parabolic equation. □ 



Prom (I33i]l . (I335]l . (l3:23]l and ([3371) . we obtain 

\\p\\m + \\pt\\m + \\u\\h3 + \\ut\\m + PWm + \\(^t\\m 

+ I {4t + ulxt + ptt + Gl + elt + + ul + 4) < c. 

JQt 

This proves Theorem 13.11 □ 



Proof of Theorem ll.lt 

Consider (|l.ip - (|1.4p with initial data replaced by (/Oq, uq, 6*9), we obtain from Theorem 
13. II that there exists a unique solution (p^, ^ 6^), such that (13.54p and (13.55P are valid when 
we replace (/>, n, 6) by (p'^, u"^, 0"^). With the estimates uniform for (5, we take 5 — )■ 0^ 
(take subsequence if necessary) to get a solution to (ll.ip - (11.4p still denoted by (p, u, 6) which 
satisfies (|3.54p by the lower semi-continuity of the norms. This proves the existence of the 
solutions as in Theorem ll.il The uniqueness of the solutions can be proved by the standard 
method like in [4], we omit it for brevity. The proof of Theorem II. H is complete. □ 



4 Proof of Theorem 11.2! 

In this section, we use the similar arguments as in Section 3 to prove Theorem 11.21 Through- 
out this section, we denote c to be a generic constant depending on pQ, uq, 6q, T and some 
other known constants but independent of 5 for any 5 £ (0, 1). 
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Denote = po + S, 9q = 9q + 6 and Pq = P{pfi, Oq), where po and 9q satisfy the 
same conditions as those in Theorem [Ol Note that ^ H^{I), Po > ^ > 0, 6*0 G H^il), 
f^x^oU=o,i = dxOo\x=o,i = 0, and 



(4.1) 



WPoWm < c, 
\\{^o)Jl-<c, 
H\\h3<c. 

Different from Section 3, we need to mollify 53. Denote = Js*g3, then G C°°{I), where 



9-i{x) = < 



-53(-x), X £ [-1,0), 
gsix), X e I, 
-53(2 -x), xG (1,2], 



and Js{-) = "^'^("^)' ^^'^ J usual mollifier such that J G C^(M), suppJ G (—1, 1), 

and /]g J{x)dx = 1. Since 53 G Hq{I), we have 53 G //q([— 1, 2]) and 



g'si-x), X G [-1,0), 
53(2;), X e I, 
g'^{2-x), xG (1,2]. 



Claim: 



gl 93 in i^H-^). as (5 ^ 0, 

\\gi\\m(l) < c||53ll/^i([-i.2]) ^ c||c/3||//i(/), for any 6 G (0,1), 
J\\/^oi93)xx\\L2(i) < c, for any 6 G (0,1). 

In fact, the proof of (j4.2p i and ()4.2p 9 can be found in [7]. We are going to prove (j4.2p q. 



(4.2) 



iy pf> - Vp^){4)xx + Vp^igi) 



3/xx 



+ Vp^igi) 



'pI + Vp 
= ^1 + A2. 

Recall Js{-) = ;^*^(;^)' conclude 

< c\\{g3)x\\L^{[-l,2]) 

< C||(ff3)^||i2(/). 

A direct calculation combining {■y/poig3)x) ^ G L'^{I) gives 

l\A2\'<c. 
By dMI), 03]) and g3]), we get (jO) .. 



(4.3) 



(4.4) 



(4.5) 
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Let Uq be the solution to the following elliptic problem for each 6 € (0, 1): 

^Oxx ~ {Po)x — Pods 7 



(4.6) 

'Ola;=0,l 



■Unl.r=n.l — 0. 



Since pf^ = + 6 e H^{I), 6^^ = 9o + 6 e -H'^(I), and e we obtain from the elliptic 

theory (see [7]), (gl]), (g^]) and (jMl) that G F^(/) n H^{I) with the following properties: 

uff no in i?^(/), as 5 0, 

< c for any 5 G (0, 1). 

Theorem 4.1 Consider the same assumptions as in Theorem Then for any T > and 
6 G (0, 1) there exists a unique global solution {p,u,6) to Ill.l\) - I^1.4^ with initial data replaced 
by (pq,Uq,6^), such that 

pGC{[0,T];H''), ptGCi[0,T];H^), pu e C{[0,T]; H^) D LH[0,T]; H^), 
Ptu^L^iQr), P>|>0, u^C{[^,T]-H^r\Hl)r\L^{[^,T]-H^), 
ut G C{[0,T];H^)nL\[0,T];H^), uu G C{[0,T]; L^) D L\[0,T]; H^), 
9eC{[0,T];H^)nL^{[0,T]-H^), 9t e C{[0,T]- H^) n L^[0,T]- H^), 
9ueL\[0,T];L^), 9>cs>0, 

where cs is a constant depending on 6, but independent of u. 

Proof of Theorem I4.lt 

Similarly to the proof of Theorem 13.11 Theorem 14.11 can be proved by some a priori 
estimates globally in time. 

For any given T G (0,+oo), let {p,u,9) be the solution to (|l.ip - (|1.4|) as in Theorem 14. 1[ 
Then we have the following estimates. 



Lemma 4.1 Under the conditions of Theorem \4.1\ it holds for any < t <T 

\\iVp)x\\L°° + \\{^/p)t\\L°° + WpWh^ + \\Pt\\H^ + WuWh'^ + llmll/fi + WVpUth^ + \\9\\h3 



+ \\VP^t\\L^ + \\p9t\\H^ + / 



{ult + pI + + Olt + P^ul + p^^) < c. 



Proof. Though the initial velocity in Theorem 14.11 (i.e. Uq) is different from that in Theorem 
13.11 (i.e. Uq), both of them are bounded in . It suffices to check if (I3.26P and (I3.39P work 
here. If do, Lemma |4. II will be obtained from p.54p . 
By ([3TT]) and gSD 

\^ut{x,Q)\ < J = - + yPoWoUox\ 



pi\9i\ + ^pi\44x\ 

l^puUO)<c. 



This gives 
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Therefore, (j3.26p is valid here. 

Multiplying (j3.7p by — — , taking t — )• 0+, and using (11. Dp *), we have 



< 



< 



+ \\f?o40ox\ + \\/?oGUcc 



+ 



+ 



+ c 



Note that \\uq^ — uox\\l'^{i) < cV5 by (ll.OP i and (14. 6p . This gives 
Therefore, (|3.39p is valid here. 

Lemma 4.2 Under the conditions of Theorem \4-l\ it holds for any <t <T 



□ 



J I Jo- 



pUtt < c. 



Proof. Multiplying ()3.23p by uu, integrating it over /, and using integration by parts. Lemma 
12.21 Lemma l4.ll and the Cauchy inequality, we have 



1 d 



2dt J I 



1 



d 



dt 



j ptUUxUt + J pttUUxUt + j I 



o3I / PtUt + TT / PttUt - -7- / PtUUxUt + / PttUUxUt + / ptUtU, 



+ / ptUUxtUt - I pUtUxUtt - j pUUxtUtt + 



d_ 
di 



Ptu 



xt 



PttUxt 



d 



1 



< ^ / ( PtUxt - -ptUt - PtUUxUt ] + C I U^t I 



+ C I Uxt + C I Ptt + - I pu^t - I PttUxt + C. 



This gives 



Uxt < ^ / {2PtUxt-ptUt 



2ptuuxut) +c J^ul^ J^pf^ + c J^'. 



'xt 



+c j pl-j^jp^-2 Pttiuxt - Pt) + c. (4.8) 
We are going to estimate the last term of the right side of (|4.8p . By (A2)-{A4), integration 
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by parts, Lemma l4.ll and the Cauchy inequality, we have 

-2 jPu{u.,t-Pt) 
= -2 j{pQ)tt{u^t - Pt) - 2 j{Pc)tt [u:,t - {pQ)t - {Pc)t] 
< -2 j [{k9x)x + ul- {puQ)^ - p9Q'ux]^iuxt - Pt) 

+c j^Pu + c jult + c j^pOl + c 

= 2 y {K.9x)t{Uxx - Px)t -'i j UxUxt{Uxt - Pt) -2 j {puQ)t{uxx - Px)t 

+2 y [pOQ'ux] , {uxt -Pt) + cj^ pI + c jult + c pel + c. 
This, combining (j3.1ip . (A2)-(j44), Lemma l4.1l and the Cauchy inequality, concludes 
-2 j^Pu{uxt-Pt) 

< c + 2y {K9x)tiput + puux)t + c j ult + c j p6l 

-2 j {puQ)t{put + puux)t + 2 j {p6Q'ux)t{uxt - Pt) + c j 

< c j^\{Kex)t? + \ j^pul + c jult + c j^pOl + c j^pl + c. (4.9) 
Substituting (jMD into (gSD, we get 

\ y PUtt + ^ y - ~[l i^PtUxt - PtUt - 2ptUUxUt) + c jult ptt + c juli 

+c j pl-jJp! + cj\{Kex)t? + cj^pel + c. (4.10) 

Integrating dilO]) over (0, t), and using ^^i, integration by parts, ([M]), dSlTT) . (1121)2, 
and Lemma l4.ll we have 




JI 



Tx)t\^ + C 



< / {2PtUxt + {pu)xUt - 2ptuUxUt) +c ul^ pft + c / / 
Ji Jo Ji Ji Jo Ji 

= / {2PtUxt - 2puutUxt - 2ptuUxUt) + c / ul^ pf^ + c / \{Hi9x 
Ji Jo Ji Ji Jo Ji 

< I [ ult + c [ p9l + c [ p^u'ul + c / pWuI + cf f ult f pl + cf [ \{Kex)t\^ + c 
^ Ji Ji Ji Ji Jo Ji Ji Jo Ji 

< I [ ^It + c [ pOl + +c [ I ul^ [ pI + c [ I \{Kex)t? + c, 
^ Ji Ji Jo Ji Ji Jo Ji 



which implies 

/ / pul + I ult<c j pOl + c j j ult j pl + cj [ \{Kex)t\^ + c. (4.11) 
Jo Ji Ji Ji Jo Ji Ji Jo Ji 

Using the Gronwall inequality and Lemma 14.11 we complete the proof of the lemma. □ 
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Corollary 4.1 Under the conditions of Theorem \4-l\ it holds for any <t <T 

Proof It follows from ([323]), Lemma O and (^2)-(^4) that 

ulxt < c / pul + c / pfuf + c Ptu^ul + c / p^u^ul 
JQt JQt JQt JQt 

2 



-C I p^U^ul^ + C \{pQ)xt? + C \{Pc)xt\ 

JQt JQt JQt 



T •'VT •'VT 

2 I ^ / ^.2 I ^ / ^2 I ^ / 11/3 l|2 j ^ / /i2 



JO JQt JQt JO JQt 

< C. 

This proves Corollary 4.1. □ 
Lemma 4.3 Under the conditions of Theorem \4--l\ it holds for any <t <T 

/ {Pxxx + Pxxt + Ptt) + / {Pxtt + '^xxxx) — 

Proof. Differentiating ()3.3ip with respect to x, we have 

Pxxxt — Pxxxx^ "^PxxxUx QPxx^xx ^^Px^xxx PUxxxx- (^•-^2) 

Multiplying (|4.12|) by 2pxxx, integrating the resulting equation over I, and using integration 
by parts and the Holder inequality, we have 

^ /" 2 _ /■ 2 f f f 

~tT I Pxxx — / Pxxx^x 12 / PxxPxxxUxx 8 / PxPxxxUxxx 2 / pPxxxUxxxx 

dt J I J I J I J I J I 

Uxx 1 1 L°° II Pxx 1 1 1 1 Pxxx 1 1 L2 
1 1 1 1 L°° 1 1 Pxxx Wl^W Uxxx II L2 ~I" 2||/5||2^oo|| Pxxx Uxxxx \ 

By Lemma 14.11 and the Cauchy inequality, we get 



J Pxxx ~ y Pxxx ^ ^ J ^^^^^ ~^ ^' i'^-^'^) 

Differentiating (j3.33p with respect to x, we have 

Uxxxx = PxxUt + 2pxUxt + PUxxt + {PxUUx)x + {pul)x + {pUUxx)x 

+ iPc)xxx + {pQ)xxx- (4.14) 

By (|4.14|) . {Aq) and Lemma HTT] we have 

^ ul^^^ <cj^ pul^t + cj^ pI^^ + c. (4.15) 
By (|4.13p , ()4.15p , Corollary 14.11 and the Gronwall inequality, we get 

pIxx < c- (4.16) 
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It follows from (fiJ5]) . (fiJ6]) and Corollary O that 

7/^ < r 

"'XXXX — ^' 

A direct calculation, combining (jl.ip i . (j3.3ip . ()4.16p . Lemma |4. 11 Corollary 14. 1 1 and Lemma 
I4.2t implies 

/ {.Plxt + +/'?*) + / pltt < C. 
Jl JQt 

The proof of Lemma 14.31 is complete. □ 

The next lemma play the most important role in getting H'^ estimates of u. 
Lemma 4.4 Under the conditions of Theorem\4--l\ it holds for any <t <T 



I 



Proof. Differentiating ()3.23p with respect to t, we have 

{pUtt)t + PttUt + PtUtt + {PtUUj: + pUtUx + pUUxt)t + Pxtt = Uxxtt- (4.17) 

Multiplying (j4.17p by p'^'^uu (72 is to be decided later), and integrating the resulting equation 
over /, we have 

- ^ j P'^^Ptuft - j [pttUt + PttUUx + 2ptUtU.jc + 2ptUUxt + 2pUtUa:t + Pxtt]{p'^^Utt) 



2 



o72+l„2_„ _ I p'ri+'^uUttU^tt-72 / P'^^ ^PxUttUxtt 



< 4iVp)t\\L^ I P^'^-^ul + C / p^^'ul- I p''Utt{pQ)xtt- I P^'Uu{Pc)xU 



+c J p'^^+^nl - J p"''^^^uuttUxtt - 272 j p^'^ '^Uxttutt{\/p)x + c 

< C ^ p^^+"2ul + C p272^2 ^ ^ p^^uMpQ)tt +12 p''~'pxUu{pQ)t 
|2 I 1 f ^72„2 ^72+2„2 ^72-l„2 



-c\\Pxtt\\h + - I P^'Ku + c I P''^^utt + c I p^'^-'utt + c 



< cj^p'^-\l + Cy^p272^2 ^ ^l^p^^ulu + cl^P^' \ipQ)u\ 
;^p272-2^2 _^ g^^|(pQ)^^|2 _^ dlp^ttlllz + C. 



Here, we have used integration by parts, the Cauchy inequality, (^2), (^3); Lemma 
Lemma |4.H Lemma 14.21 and Lemma 14.31 

After the third term of the right side is absorbed by the left, we have 

d 



dt 



j p'<-+'ul + / p'<-ul,, < c [ p^--^ul + c / p^^-ul + c\ p^- \^pQ\t? 
Jl Jl Jl Jl Jl 

+cj^p^^^-^ul + cJ^p\ipQ)u\^ + c\\p,u\\l2 + c. (4.18) 
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By Lemma 14.21 we know Jq^ pUf^ < c. This implies that the minimum of 72 we should take 
in (j4.18p is 2. Substituting 72 = 2 into (|4.18p . we have 



di 



■ / p\l + / p^lu <c [ pul + c [ p\ ipQ)uf + c [ pIu + c. (4.19) 
Ji Ji Ji Ji Ji 

We are going to estimate Jj p\{pQ)tt\'^ ■ Using Lemma \^A] (A4) and Lemma HTSl we have 

J^p\{pQ)ttf = j^p\pttQ + 2ptQ'et + pQ"el + pQ'eu\^ 

< c + c\\et\\lo. + c I p^el. (4.20) 



I 



Substituting ()4.20p into ()4.19p . integrating the resulting inequality over (0,t), and using 
Lemma |4.H Lemma 14.21 and Lemma 14.31 we get 

/ P^ul + [ [ p^ulu < [ p\l{x,0) + c. (4.21) 
Ji Jo Ji Ji 

Using (gl]), gZD, M, (I3:23]l and g^Ds, we have 

P^ulix,0) < c, 

which combining (|4.2ip completes the proof. □ 
Lemma 4.5 Under the conditions of Theorem \^.l\ it holds for any <t <T 

/ PuLt + / {uLxt + P^lxt) < c- 
J I JQt 

Proof By (j3.23p . we have 



pul^t < c p^uft + c / ppfuj + c / ppW^I + / P^u^ul 



I 



+C p'^U^ul^ +cj^p\ {pQ)xt\^ + / 



2 



where we have used {A2)-{Ai), Lemma |2.2|, Lemma |4.1|, Lemma 14.21 and Lemma 
It follows from (j3.37p and (^5) that 



+ c/ p\K'\Hlel 



[ pel,, < c[ p\K'\^9^el, + c [ p\K"\^9^et + c [ 

JQt JQt JQt JQ'. 

+c [ pWel + c [ p'\Q"\X + cj ppI\q'\H', 

Jqt Jqt Jqt 

+c p\ipuQ'9x)t\'^ + c p\{pOQ'ux)t\'^ + c pulult, 
Jqt Jqt Jqt 



which, combining (^4), (^5) and Lemma l4.ll gives 

cT 



[ P^lxt < c[ 9lt + c[ p^9l + c[ \\9t\\l^+ [ ult 

JQt JQt JQt Jo JQt 



< c 9lt + c i 


f P'el + cf 


JQt J 


Qt Jo 


< c. 





(4.22) 
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Differentiating (jS.SSp with respect to t, we get 

Uxxxt = '^{■\/p)xy/pUtt + pUxtt + PxtUt + PtUxt + PxtUUx + Ptul + PtUUxx + PxUfUx 
+2pUxtUx + pUtUxx + PxUUxt + PUUxxt + {pQ)xxt + {Pc)xxt- 

This, together with (j4.22p . (Aq), Lemma 12.21 Lemma l4.ll Lemma 14.21 Lemma 14.31 Lemma 
14.41 and Corollary 14.11 implies 

/ ulxxt <c + c [ {pel^t + elt) + c [ ||6't iiioo + c / < c. 

JQt JQt Jo Jqt 

This completes the proof. □ 

Prom (j4.15p . Lemma 14.31 and Lemma 14.51 we get the following corollary immediately. 
Corollary 4.2 Under the conditions of Theorem \4-l\ it holds for any <t <T 

< r 

"'xxxx — ^' 

Corollary 4.3 Under the conditions of Theorem\4-l\ it holds 



92 < c 

xxxx — 



'Qt 

Proof. Differentiating (j3.42p with respect to x, we have 

i^Oxxxx = -4:K'ex9xxx - 3k - 3k"919xx - {n"9x)^ - 2 {uxUxx)x + [pQ'Gxt).^ 
+ {pxQ'Ot), + {pQ"0xOt), + ipuQ'9x)xx + ipOQ' Ux^xx • 

This, combining (^5), (Aq), Lemma l3.3t Lemma l4.ll and Lemma l4.5t implies 



^Ixxx <c + c I {p6l,, + O + c [ + c I \k"'\^91 < c. 

Qt JQt Jo JQt 

This proves Corollary 4.3. □ 

Lemma 4.6 Under the conditions of Theorem\4- 1\ it holds for any <t <T 



Pxxxx ~^ I ^xxxxx — 
JQt 

Proof. Differentiating (|4.12p with respect to x, multiplying the resulting equation by Ipxxxx-, 
integrating over /, and using integration by parts, Lemma [4.H Lemma [4.3|, Corollarv 14.21 and 
the Cauchy inequality, we get 



d 
di 



2 _ f 2 f f 

Pxxxx — ^ I Pxxxx^x 20 / Pxxx PxxxxUxx 20 / PxxPxxxx'^xxx 

Jl Jl Jl 

10 / PxPxxxxUxxxx 2 / pPxxxxl^xxxxx 

Jl Jl 

< c ^ pIxxx + c ^ ulxxxx + c. (4.23) 
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Now we estimate the second term of the right-hand side of (j4.23p . 
Differentiating (j4.14p with respect to x, we have 



PxxxUt + SpxxUxt + ^PxUxxt + PUxxxt + {PxUUx)xx + (pUx) 



XX 



~\~{pUUxx)xx xxxx xxxx ■ 

This, combining {Aq), Lemma 12.21 Lemma l4.ll Lemma 14.31 and Corollary 14.21 concludes 

/ ul^^^^ <c I ul^t + c I ul^^t + c pI^^^ + c el^^^ + c. (4.24) 
Ji Ji Ji Ji Ji 

Substituting (j4.24p into (|4.23p . and using Corollary 14. 1[ Corollary 14.31 Lemma and the 
Gronwall inequality, we obtain 

pIxxx < c. (4.25) 



It follows from l![^7Z^ , ([O^ , Corollary HTTl Corollary [O] and Lemma [13] that 



'^xxxxx ^ C. 



This completes the proof of Lemma 4.6. □ 
Corollary 4.4 Under the conditions of Theorem \4-i\ it holds for any <t <T 

/ {pltt + pixxt) + / {Pm + plxtt) < c. 
JI JQt 

Here we haye used the following inequality when we get the upper bound of put- 

pWtt = 2 [{y/p)x^/p\^ ul < cpul. 
From the above estimates, we get 

\\{Vp)x\\l°° + \\{Vp)t\\L'=° + \\p\\m + \\pt\\H» + WpttWm + 

3 

+11^411^:^1 + Wp^uuWl^ + W^Uxxth^ + PWm + WVpGtWi'^ + Wp^xtWi^ 

1 / 2 2 2 2 2 2\ 

+ / [P "^xtt + P'^tt + '^xxt + '^xxxt + '^xxxxx) 
JQt 

+ / {pit + plxtt + Olxxx + + Olt + pelxt + < c. (4.26) 

JQt 

From (14:26]) and (|335]) . we get 

WpWh* + WPtWn'^ + WpttWm + WuWh^ + \\ut\\H-^ + \\utt\\L^ + WOWn'^i + Pt\\m 

+ / (""xtt + '^xxxt + ''^xxxxx) + / [Pttt + Pxxtt + ^xxxx + ^xxt ~^ ^tt) ^ c((5), 
JQt JQt 

where c{5) is a positive constant, and may depend on 5. 

The proof of Theorem 14.11 is complete. □ 
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Consider ()l.ip - (jl.4p with initial data replaced by (pg, Ug, 9q), we obtain from Theorem 
14. II that there exists a unique solution {p\ u\ e^) such that (14261) and ([335]) are valid when 
we replace (p, u, 9) by {p^ , u^, 9^). With this estimates uniform for 6, we take (5 — t- 0^ ( take 
subsequence if necessary) to get a solution to (|l.ip - (|1.4p still denoted by (p, u, 9). By the 
lower semi-continuity of the norms, we have 

\\{^/p)x\\l°° + ll(\/p)t||L°° + \\p\\m + WpiWh^ + WpttWm + + \\ut\\m 

+ \\^/pU^xt\\L^ + \\9\\m + \Wp9t\\L^ + \\p9xt\\L^ + / {ulxt + ^Ixxt + ^Ixxxx) 

JQt 

//'2 2 2 22\ 

+ / [Pttt + Pxxtt + (^xxxx + +(^xt) < C, 

which proves the existence of the solutions as in Theorem [L2j The uniqueness of the solutions 
can be proved by the standard method like in [3], we omit it for brevity. The proof of Theorem 
11.21 is complete. □ 
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A blow-up criterion of strong solution to a 3D viscous 
liquid-gas two-phase flow model with vacuum 
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In this paper, we get a unique local strong solution to a 3D viscous liquid-gas two-phase 
flow model in a smooth bounded domain. Besides, a blow-up criterion of the strong solution 
for yyu > /I is obtained. The method can be applied to study a blow-up criterion of the strong 
solution to Navier-Stokes equations for ^jj. > A, which improves the corresponding result about 
Navier-Stokes equations in Iil5j where 7/i > A.. Moreover, all the results permit the appearance of 
vacuum. 

Keyword: Liquid-gas two-phase flow model, local strong solution, blow-up criterion, vac- 
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1 Introduction 

In this paper, we consider the following 3D viscous liquid-gas two-phase flow model 



Huanyao Wen^'^* Lei Yao^]^ Changjiang Zhu'* 



Abstract 



uum. 




{mu)t + div(mM ^u) + VP{m, n) - /jAu + (p. + A)Vdivu, in Q X (0, oo), 



(1.1) 



with the initial and boundary conditions 



{m,n,u)\t=o ^ {mo,no,uo), in O., 



(1.2) 



u{x, t) = 0, on dD. X [0, oo), (1.3) 



where Q c is a smooth bounded domain. Here m = aipi and n - OgPg denote the hquid mass and 
gas mass, respectively; /u, A are viscosity constants, satisfying 



yu > 0, 2iJ + 3A>0, 



(1.4) 
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which impUes fi + A > ^fx > 0. 

The unknown variables at, Og € [0, 1] denote respectively the liquid and gas volume fractions, 
satisfying the fundamental relation: ai + ag = 1. Furthermore, the other unknown variables p/ and 
Pg denote respectively the liquid and gas densities, satisfying equations of state: pi - pi^ + ^ C'" , 

a, 

Pg = ^, where ai, Ug are sonic speeds, respectively, in the liquid and gas, and Pifl and pio are the 
reference pressure and density given as constants; u denotes velocity of the liquid and gas; P is the 
common pressure for both phases, which satisfies 

P(m, n) - C° {-b{m, n) + ^Jb{m, n)^ + c{n)) , (1.5) 
with C° = ia2, ^0 = P/,0 -^>0,ao = C^f and 

b{m, n) = k()- m-\ — n = ko - m - aon, 
\a,j 

lag\^ 

c(n) - 4ko — n - 4kQaQn. 
\a,j 

For more information about the above models, please refer to lfm[T4ll20ll and references therein. 

The investigation of model (11.11 ) has been a topic during the last decade. There are many results 
about the numerical properties of this model or related model. However, there are few results pro- 
viding insight into existence, uniqueness, regularity, asymptotic behavior and decay rate estimates 
concerning the two-phase liquid-gas models of the form (11.11) . Let us review some previous works 
about the viscous liquid-gas two-phase flow model. For the model (11.11 ) in ID, when the liquid is 

incompressible and the gas is polytropic, i.e., P(m, n) = Cp] I — ^ — ) , Evje and Karlsen in Q stud- 

\pt - m) 

ied the existence and uniqueness of the global weak solution to the free boundary value problem with 

p - p{m) = k\ —, jS € (0, t), when the fluids connected to vacuum state discontinuously. 

{pi-rnf^^ 

Yao and Zhu extended the results in |4 j to the case fi e (0, 1], and also obtained the asymptotic be- 
havior and regularity of the solution, see lITSl . Evje, Flatten and Friis in 121 also studied the model 

with p = p{m, n) = k2 ^ (J3 e (0, ?)) in a free boundary setting when the fluids connected to 

(pi - nif^^ 

vacuum state continuously, and obtained the global existence of the weak solution. Also, for the case 
of connecting to vacuum state continuously, Yao and Zhu investigated the free boundary problem to 
the model with constant viscosity coefficient, and obtained the existence and uniqueness of the global 
weak solution by the line method, where a new technique was introduced to get the key upper and 
lower bounds of gas and liquid masses n and m, cf. |[T9l . Specifically, when both of the two fluids 
are compressible, one can consult the reference ||3]. For multidimensional case, the results are few. 
Recently, Yao, Zhang and Zhu obtained the existence of the global weak solution to the 2D model 
when the initial energy is small, see L20 I . Furthermore, they proved a blow-up criterion in terms of the 
upper bound of the liquid mass for the strong solution to the 2D model in a smooth bounded domain, 
cf. 1211 . Because of the complexity of the pressure P{m, n), they in 11211 can only deal with the case: 
there is no initial vacuum, i.e., wiq > 0, «o > 0. Then, what will happen when the vacuum appears? 
In this paper, we prove the local existence of strong solution and give a blow-up criterion to the 3D 
viscous liquid-gas two-phase flow model in a smooth bounded domain with vacuum. 

The main results are stated as follows. 
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Theorem 1.1 (Local existence). Let D. be a bounded smooth domain in and q € (3, 6]. Assume 
that the initial data mo, «o. satisfy mo, «o ^ W^'^(n), mo ^ Hq{Q.) n H^(Q.), < ^gmo < «o ^ ^ow^o 
Q, where £q andlo are positive constants. The following compatible condition is also valid: 

- fiAuQ - (ju + /l)VdivMo + ^P{mo, «o) - V"^o5' ■^o'we g € L^(Q). (1.6) 

r/je«, ?/jere exist aTo > and a unique strong solution (m, n, u) to the problem di.iD -f TOl) . 5mc/z that 

0<Sj^m<n< Tom, (m, «) e C([0, Tq]; lV''^(f2)), (m„ n,) e L"(0, To; L^(Q)), 

P e L~(0, To; lV'''''(n)), u e L"(0, To; //|J(Q) n //^(Q)) n 1^(0, To; W'^-'^ia)), 

^J^Ut € L~(0, To; L^m, Ut e 1^(0, To; //^(Q)). (1.7) 

Furthermore, under the assumption 

25 

A < —fi, (1.8) 
we can establish a blow-up criterion of the strong solution: 

Theorem 1.2. Under the assumptions of Theorem \l.l\ ifT* < oo is the maximal existence time for 
the strong solution (m, n, u){x, t) to the problem di-iP -f TOl) stated in Theorem \l.l\ then 

lim sup \\m\\L->(oj-L'^(ii)) = oo, (1.9) 

provided that di.^l ) holds. 

Remark 1.1. (i) For Q. = R^, we can also get a unique strong solution to (ll.ll )- (ll.5l ) and the blow-up 
criterion (ll.9l l by using the ideas of lH] [I5l to modify the proofs of Theorem 11.11 and Theorem 11.21 
shghtly. 

(ii) The proof of Theorem 11.11 and Theorem 11.21 implies that the following blow-up criterion 
would be obtained if the restriction (ll.Sl l is removed: 

hm sup (||m||z,~(o,r;L~(n)) + II ^uWj^HOTu'm) = °°' (l-^O) 

T->T' 

where ^ + ^ < I and 3 < < oo. (11.101) is similar to Q. 

(iii) Under the assumption (11.81 ). we can use our methods in Lemma \5?2] together with the esti- 
mates in 1 15 1 to get the following blow-up criterion of strong solution to Navier-Stokes equations: 

lim sup \\p{t)\\L'-{0,T;L'-m = «3- 

This relaxes the restriction l/u > Ain (151 . And our result can be viewed to be a generalization of |[T5l 

We should mention that the methods introduced by Sun, Wang and Zhang in (15), Cho, Choe and 
Kim in [IJ for Navier-Stokes equations will play a crucial role in our proof here. There are many 
results about blow-up criterion of the strong solution for the Navier-Stokes equations in addition to 
Q. For the 2D compressible Navier-Stokes equations. Sun and Zhang in [16] obtained a blow-up 
criterion in terms of the upper bound of the density for the strong solution. For the 3D compressible 
Navier-Stokes equations. Sun, Wang and Zhang in [15 ] obtained a blow-up criterion in terms of the 
upper bound of the density for the strong solution, under the restriction A < Ip. In both papers, the 
initial vacuum (po > 0) was allowed and the domain included both the bounded smooth domain and 
R^, A'^ = 2, 3. It also worths mentioning recent works ||8]|9], under the assumptions 

N = 2, p>0, p + A>0, a = T^\ 
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or 

N - 3, A < 7/1, // > 0, and 2/i + 3/1 > 0, Q is a smooth domain including R , 

Huang and Xin proved the following blow-up criterion: if T* < oo is the maximal time of the existence 
of the strong solution, then 

lim r \\VumL-(n)dt = oo. (1.11) 
T-*'' Jo 

Huang, Li and Xin in their recent paper [10 | removed the restriction A < l/u for N = 3, and got the 
blow-up criterion of strong solution: 

lim r ||D(M)(OllL~(n)^/? - oo, 

T->T' Jo 

where D{u) - ^(Vm -i- VmO- For the non-isentropic compressible Navier-Stokes equations, under the 
conditions: N = 2, /u > 0, ju + A > 0, Q = or [0, 1]^; N = 3, A < Ifi, /u > 0, and 2iu + 3A>0, O 
is a smooth bounded domain, please refer to |[T2l l5l. 

In Theorem 11.21 we give a blow-up criterion in terms of the upper bound of the liquid mass 
under the relaxed restriction (II. St . which improves the corresponding result about Navier-Stokes 
equations in [15] where 7yu > A. Here, if the liquid mass is upper bounded, we can obtain a high 
integrability of the velocity, sup j m\u\''dx < C, for some r € (3,4], see Lemma [S!2l Moreover, in 

Q<t<T 

order to overcome the singularity brought by the pressure P{m, n) when there is vacuum, we need the 
assumption: < Sj^niQ < «o < ^o'mo, where Sj^ and Iq are positive constants. 

2 Preliminaries 

In this section, we give some useful lemmas which will be used in the next three sections, where 
A^ = 2,3. 

Lemma 2.1. Let Q. c R'^ be an arbitray bounded domain with piecewise smooth boundaries. Then 
the following inequality is valid for every function u e W^'^iQ) or u e 

ll«llL/(n)^Ci||V«||2,(„)||«||[;«^^, (2.1) 

where a = (1/r' - l/p'){\/r' - l/p + l/N)'^; moreover, if p < N, then p' e [r',pN/{N - p)]for 
r' < pN/{N - p), and p' e [pN/iN - p),r']for r' > pN/(N - p). If p > N, then p' € [r',oo) is 
arbitrary; moreover, if p > N, then inequality ([2.1\l is also valid for p' = oo. The positive constant Ci 
in inequality f l2.iD depends on N, p, r', a and the domain Q. but independent of the function u. 

Lemma 2.2. Let Q. c be an arbitrary bounded domain with piecewise smooth boundaries. Then 
the following inequality is valid for every function u € W^'P{Q.): 

\M\u"in) ^ C2(||M|L.(fi) + \m\l,^a)M[''W' (2.2) 

where N, p, r', p' and a are the same as those in Lemma 12.71 The positive constant C2 in inequality 
l \2.21) depends on N, p, /, a and the domain Q. but independent of the function u. 

The above two lemmas can be found in |[T3l [171 and the references therein. 
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Next, we give some V\p e (1, oo)) regularity estimates for the solution of the following boundary 
problem: 



LU ■- fiAU + + A)VdivU = F, in Q, 

U{x) - 0, on da. ^ ' 

Here Q c is a bounded smooth domain, L is the Lame operator, U = {U\, U2, • • • ,Un),F = 
{F\,F2, • ■ • ,Fi\i). From (ll.41 i. we know that (12.31 ) is a strong elliptic system. If F e then 
there exists an unique weak solution U € //qC^I). In the subsequent context, we will use L^^F to 
denote the unique solution U of the system (12.31 ) with F belonging to some suitable space such as 
W~^'P{a.). Sun, Wang and Zhang in |[T5l[T6l give the following estimates: 

Lemma 2.3. Let p e {1, 00), and U be a solution of d2.3| ). Then there exists a constant C depending 
only on ^, A, p, N and Q. such that 

(1) ifF e LP{Q.), then 

\\U\\w2.n(a) < CWFhrisiy, (2.4) 

(2) ifF € W-''P(n) (i.e., F = div/ with f = {fij)NxN, fij e LP{a)), then 

l|f/|lwi./'(n) < C||/||L.(n); (2.5) 

(3) ifF = div/ with fj = dkh), and h), € W,l'^(a)/or i, j,k ^ 1,2,- ■■ ,N, then 

J Ij IJ U 

\\U\\LP(n) < C\\h\\Lnn). (2.6) 

Lemma 2.4. If F ^ div/ with f - {fij)NxN, fj e ^"(H) n L^{Q.), then Vf/ e BMOiQ) and there 
exists a constant C depending only on p, A and Q such that 

\\^U\\BMO(m ^ C(||/||L~(n) + ll/llL2(n))- (2.7) 

Here BMO(Q.) denotes the John-Nirenberg's space of bounded mean oscillation whose norm is de- 
fined by 

ll/llBMO(n) = ll/llL2(n) + [f]BMO{Ci), 

with the semi-norm 

[f]BMO(Si) = sup f \fiy) - fiiAx)\dy, 

xeCl,re(0,d) Jn^x) 

where Q.r(x) = Br{x) n Q, Br{x) is the ball with center x and radius r and d is the diameter of fi. For 
a measurable subset E of R^, l^l denotes its Lebesgue measure and 

/n,(x) = r f(y)dy = -r^ I f(y)dy. 

Lemma 2.5. Let Q. be a bounded Lipschitz domain in R^ and f e W^'P(Q.) with p € (A'^, 00). Then 
there exists a constant C depending on p, N and the Lipschitz property of the domain Q. such that 

ll/llL"(n) < C (1 + WfWBMOiSi) ln(e + ||V/|b.(n))) . (2.8) 
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3 Global existence for the linearized system 

Consider 

fUt + div(mv) = 0, 
Ht + div(«v) - 0, 

mut + mv ■ Vm + VP(m, n) - Lu, in Q. X (0, oo), 
with the initial and boundary conditions 

(m,n,u)\t=o ^ {mo,no,uo), in Q, 



(3.1) 

(3.2) 
(3.3) 



u{x, t) = 0, on dQ. X [0, oo), 

where Q c is a smooth bounded domain. 

Throughout the rest of the paper, we denote W'^''^ - W'^'P{Q) for > and I < p < oo with the 
norm || • \yk.p. Particularly, H'' - W*"'^, and LP - W'^'P. Qt ^Tix [0, T]. 

Theorem 3.1. Let D.be a smooth bounded domain in B? and q e (3, 6]. Assume mo, n^ € W^'^, uq € 
//I n H^, mo >6>0,no>6>0,ve C([0, T]; //^ n H^) n L^{0, T; W^'?) and Vt € L^{0, T; Hi). Then 
there exists a unique strong solution (m, n, u) to dj.iP -f TOl) such that 



(m, n) e C([0, T]; W^''?), (m„ n,) e C([0, T]; L'?), 

P € C([0, r]; W^-''), M e C([0, T]; //^J n Z/^) n L^{Q, T; W^-'?), 

M, e C([0, rj; L^) n 1^(0, T; H^), m>0,n>Oin Qt. 

Proof. By Q, (0)1 and (I3D2, we get 

'm,n e C([0, T]; W'-'?); m„n; e C([0, T]; L'?), 

sup IWOIIwi" < l|mollwi'!exp|c \\Vv{s)\yi.<,ds\ , 
0<t<T ^ I 

sup IKOIIw'" ^ INIIw'^expIc f, ||Vv(s)||iyi.,<i5l, 
o<t<T ( I 

0<6exp {- ||Vv(s)||l» Js} < m < ||mo||L~ exp { ||Vv(s)||l» Js} , 
< <5exp{- \\Vv{s)\\L-ds} <n< ||«oIIl~ exp jX'' ||Vv(^)||i» Js} . 

This immediately gives 

P{m,n) e C([0, T]; W'-''), /'(m,«); e C([0, T]; L'?). 
It follows from (ESs, dMll, and Q that 

M e C([0, r];//^ n //2) n L^{0, T; W^'^), € C([0, r];L2) n 1^(0, T;//^). 



(3.4) 



(3.5) 



4 Proof of Theorem [LI 

In this section, we get a unique local strong solution to (ll.ll )- (ll.5l ) with mo > 6 > 0, no > 5 > 0, 
and obtain some estimates uniformly for 6 (see Theorem 14. 11 1. Theorem 1 1.1 1 will be obtained after 
constructing a sequence of approximate solutions (m*^, n'', u^) by giving the initial data (mo, no) 
in Theorem 11.11 a lower bound 6, using the estimates in Theorem 14.11 and taking 5^0 (taking 
subsequence if necessary). 
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Theorem 4.1. Under the conditions of Theorem U . 1\ we assume mq > 6 > 0, no > 5 > 0. Then there 
exists a time Tq > independent of 5 and a unique strong solution (m, n, u) to di.iD -f TOI) such that 

(m, n) e C([0, To]; (m^, e C([0, Tq]; L«), 

P € C([0, To]; W^-'?), M e C([0, Jq]; //^ n Z/^) n 1^(0, To; W^'''), 
e C([0, ro];L2) n 1^(0, Tq;//^), m>0, n>Oin Qt^. 

Moreover, we have the following estimates: 

{ sup LmN2 + /;°/^|V«,p<C, 

0<t<To 

ll"llL2(0,ro;w2.'?) + ll"llL»(o,ro;//in//2) + \\'^h"'{0,To;W^-i) + ll«llL~(o,ro;Wi-9) ^ C", 
^ - lo"^ ^ n < som, in Qtq, 
\\P\\l'"(Qtq) + ll^mllL°°(ero) + ll^;illL"(eTo) ^ C, 
where C is a positive constant, independent of 6. 

To prove this theorem, we first construct a sequence of approximate solutions inductively as 
follows (similar to yj): 

(i) Define m° = 0, and assume u^-^ € C([0, T]; //^ n Z/^) q /^2(-q^ 7.. ^2,<?) p ^i^g^ //^) ^as 
defined for A; > 1. 

(ii) By Theorem l3.1[ we can get {m^, n^, u^) with the regularities in Theorem 13. 1 I satisfying 

+ divCm^^M*'"') = 0, 

n\ + di\{n^u^-^) = 0, (4.1) 
m'^wf + m^u^-^ ■ Vu^ + VP^ = Lu'', in Q X (0, T], 

where = P{m^, n^). The initial and boundary conditions are stated as follows 

{m'',n^,u'')\t=o = {m,m,uo), in Q, (4.2) 



/(x, t) - 0, on dO. X [0, r]. (4.3) 



Throughout this paper, we denote 



0^(0 = max (1 + \\m\t)\\LA , *F^,.(0 = max (l + f 
i<A</f^ > \<k<K\ Jn 



for r € (3, 4] and K e Z+. The next step is to make some estimates for (m*^, n'^, u'^) (k > 1) independent 
of k and 6. 

Lemma 4.1. Under the conditions of Theorem \4.1\ we have for all k > 1 

< Som*" <n''< lom'^, in Qt- 
Proof It follows from H.lh and (14.11 )9 that 

nM k^i -In'' 



+ • V — =0. 



/, \ m^ 



This implies 



-f-f4l(^(^;^'0,^) = 0, (4.4) 



ds \m' 
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where X{s; x, t) is given by: 

—X{s; X, t) = M*"' {X{s; x,t),s), < s < t, 
as 

X{t; X, t) - X. 

Integrating (14.41 ) over (0, t), and using the assumption s^mq < uq < lomQ, we complete the proof of 
Lemma 4.1. □ 



Lemma 4.2. Under the conditions ofTheorem \4.1\ we have for all I < k < K 

0<P'' < CO/f(0 in Qt, (4.5) 

< < C in Qt, (4.6) 

< P^, < C in Qt, (4.7) 
where C is a positive constant, independent ofK, 6 and T. 
Proof. (14.51 1 can be obtained by Lemma |4~T] and (11.51 ). A direct calculation gives 



P^, ^ C° -! 1 - ' \ > 0, (4.8) 

[ ^Jb^{m'',n^) + c{n'') ] 

P\, - C" j ao + ={m'' + a^n^ + ^o) i > 0. (4.9) 

Obviously, we get (14.61 ) by (14.81) . To get (I4.7I ). it suffices to prove 



b{m^,n^) 



k k 1 \ 

m + «()« + Kq 
, ^b^{m^, «*) + c{n^)j 



< C. 



In fact, 

+ a^n'' + Icq ]^ (m'^f + al{n^f + kl + Ikom'' + lam^n^ + laok^n^ 



, yjb'^{m^,n^) + c{n'') ) {m^)^ + a^in^)^ + ~ "^^om^ + la^m'^n'' + laok^n'' 

AkQWi'' 



= 1 + 



< 1 + 



(ko - m*)^ + a^in^^)^ + lagm'^n'^ + la^k^n^ 
4kom'' 



la^k^n^ 
< C, 

where we have used Lemma |4TT] This completes the proof of Lemma 4.2. 
As in lfT5l . we denote w'^ - u!^ - h'^, where /j^ is the unique solution to 

ilh'' = VP'', innx{0,T], 
\h%a = 0. 

From Lemma 1231 we get for any p € {l,oo) 

i\\h'\\wu> < CWP'Wlp, 
\\\h'\\w2.,> < C\\VP'\y. 
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(4.10) 



(4.11) 



SJ)3, €3 and KM imply 



jm^wf - Lw'' = m*^F^ in Q X (0, T], 

w\x, 0) = uo- L-^VP{mo, no), in Q, (4.12) 
\w%ci = 0, 



where 



Lemma 4.3. Under the conditions of Theorem \4.1\ we have for all I < k < K, 3 < r < 4 
Jo Jn Jn I Jo J 

r r IV VP <C sup (D^COexp |c r [Ojf(^)]^['Fi,,,-(^)]^^/4' 
Jo Jn o<?<r I Jo J 

where C is a positive constant, independent ofK, 6 and T. 

Proof. Multiplying (14.121 ) by wf , integrating over Q, and using integration by parts and Cauchy in- 
equality, we have 

Jn ^ Jn 

< \ [ m'lw^f + i r mW, 
^ Jn ^ Jn 

which imphes 

r m'=|wfp + -^ r [;i|Vw*^P + (/I + i)|divw'^p] < f m*^|F*^p. (4.13) 
Jn Jn Jn 

Now we estimate the term of the right side in (|4.13l l as follows: 

fm'\F'\^ 
Jn 

< C r m*^|M*^-ip|V/p + C(l)if(0 r |L-^Vdiv(PV-i)p + C(l)if(0 f (|m*^p + |P*^p)|divM*^"'P 
Jn Jn Jn 



< C 



< C 



-'n 

r-2 



Jn J [Jn J Jn Jn 

2 

Jn J Jn Jn 



(iivw^ip 2. + m'f ) + c[<Djf(o]' r [ 

\ L7^ LF^/ Jn 



< C[(i>K{t)]^U¥KAt)]' 2r + IIV/jIP 2r 1 + C[<Djf(0]' I [//IVm^^-^P + (ju + i)|div/-ip] 



+C[OKit)]^ f [p\Vu''-^f + (p + A)\diyu''-^f], (4.14) 
Jn 



where we have used Lemma l2!2l Lemma 1231 Lemma \47\] Lemma |42] and Young inequaUty: ab < 
saP + (sp) p q b'^ for any £ > 0, p, ^ > and ^ + ^ - 1- 
By Lemma 123] and (14.121) . we have 

IIV^w*^!!^, < CO^(0 r m^lwf |2 + CO^(0 r m'\F'\\ (4.15) 
Jn Jn 

Substituting ( 14.151) into ( I4.14D . we have 

< CmKiOf-^l'VKAmH f m%\\^+ r m'^lf*^^] 

+C[O^(0]'^[^if,.(0]^(£^ + i) r ti\^w''\^ + C[^K{t)f^Vi'KAt)V 

Jn 

+C[<D^(0]^ r MVw'^-'P + Ou + ^Mvw^^-'Pl+CLO^CO]^ r MV/j^-ip + (/i + ^)|div/j*^-ift 
Jn Jn 



\Jn Jn 



n 

ka\ 



+C[OK0]'^[^i?,.(0]^(e^ + l) r //|Vw^P + C[O^(0]^[Yi,,,-(0]' 

Jn 

+C[O^(0]' r bulVw'^-'P + + ^)|divw*^-'p] + C[O^(0]^ 
Jn 



;n 

where we have used (14. lit . 
X 

AC 



Take £ - ^ [^Km ^ Vi-K^m ^ , we have 



Jn J Jn Jn 

+C{^K{t)f \ |ju|Vw^~'P + (/i + i)|divw*^-ip]. (4.16) 
Jn 

Combining (14.131 ) and (14.16b . we get 

^ r m*^|wfP + 4 r l;/|Vw*^P + Ou + /l)|divw*^p] 
3 Jn Jn 

Jn 



+C{^K{t)Y [//IVw^-^P + (// + i)|divw'=-T]■ 
in 



/i 

Jn 

Integrating over (0, t), we have 

nm*|wfp+ r |;i|Vw*^P + (ju + ^)|divw^p] 
! Jn 

< c r [a)if(^)]5^[^^,.(^)]A r ^iivw'^p + c r [a)^(^)]'['F^,,-(^)]^ 

Jo Jn Jo 

+C r [<l)7f(s)]^ r [/i|Vw'^"^P + (yU + /l)|divw'^"'p] + C. (4.17) 
Jo Jn 
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Denote Axit) = sup f [ju\Ww \^ + (p + /i)|divw p], we obtain from (14.171 ) and noticing ^ > 3 for 

l<k<K 

3 < r < 4 and (D/f > 1, ^'Kr > 1 



By Gronwall inequality, we get 



Jo 



Jo 



Jo 



+C [<Di,(T)]^[^^,,(T)]-exp<^C [<^K(.s)]-m'K,M]-ds}dT 



(4.18) 



where we have used the inequality: y < exp{j} for j > 0. By (14.151 ). (14.161) . (14.171) and (14.181) . we 
complete the proof of Lemma 4.3. □ 

From (14.111 ) and Lemma 1431 we immediately give the following corollary. 

Corollary 4.1. Under the conditions ofTheorem \4.1\ we have for all \ < k < K and 3 < r < 4 

l|V/|lL»(or L2) <C sup Oj,(0exp |c f [0^(s)]^[^^,.(s)] A , 

Q<t<T { Jo ) 



and 



sup O/f(0 + Vf 

0<t<T 



f 

Jo 



exp^C I [OK{s)]^U¥KAs)]^ds\, 



l|VM''llL2(0,r;L6) < C 

where C is a positive constant, independent ofK, 5 and T. 

Now we give higher order estimates for m*. 
Lemma 4.4. Under the conditions ofTheorem \4.I\ we have for all \ <k < K and 3 < r < 4 



r m^iM*^p+ r r (//ivM*^p+(//+i)idivii*^p) 

Jn Jo Jn 

Jo 



< C sup {^K{s)f exp C [<Dif (5)] - Vi'KAs)]-ds + 
o<i<r I Jo J -^n 



Jn 



+C sup ^k{s) r <1)^(t)(i+ V^m* ,)||Vm*^||2 exp |c f ^^/^ 

0<4<T Jo ^ I Jo 

+C ^up^^K{s) ^k{t)(i + II V^M*^"^||^,)l|V/-i||2,exp|c^ [0^(^)]^['l'jf,rW]^^/s 



where iJ^ - + ' • Vm*^, < t <T, and C is a positive constant, independent of K, 6 and T. 
Proof. (14. Ih can be rewritten as 

m'^ii'' + VP'' - Lu'' = 0. 
Differentiating (14.191 ) with respect to t, and using (I4.1l )i . we conclude 

m'^u'l + m'^u'"'^ ■ lit + IP] + divCVP'^ ^u^^^) 



(4.19) 



- Lu^ - L(it~^ ■ Vm*^) + div(LM'= ®u^~'). 



k ^ ,.k-\-, 



(4.20) 
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Multiplying (14.201 ) by ii^, integrating the resulting equation over Q, and using integration by parts, we 
obtain 

= r (Pfdivii*^ + {u^-^ ■ Vm*^) • VP*) +^i \ V(m*^"' • Vm*^) : Vm*^ + (// + i) f div(M*^-' • VM*^)divM*^ 
Jn Jo. Jn 

- (m*-^ • Vm*^) • (/iAm*^ + (yu + /l)VdivM*^) 

= J (Pfdivii*^ + (/-^ • Vm*^) • VP*) +^l ^ [V(/-i • V/) : Vm* - (m*"' • Vm*) • A/] 

+(// + /{) ( [div(M*-^ • VM*)divM* - (m*-^ • Vm*) • Vdivw*] 
= h+h + h- (4.21) 
Now we estimate I\, h and /a as follows: 

h = £[(P^,mf + P>f)divM* + (M*-i-VM*)-VP*]jx 

= J [i-m'^Pl, - ?i*P*JdivM*"MivM* - m*"^ • VP^divw* + (m*"' • Vii*) • VP*] 

= J [i-m'^P'lt - «*P*JdivM*-MivM* + P*div(M*"MivM*) - P*div(M*"' • Vii*)] 

= J [(-'w^'Pi. - «*^P^JdivM*-idivii* + P* (divM*-Mivii* - (V/'^)' : Vii*)] 
< Ccl)^(OI|V/-i||i2||VM*||i2, (4.22) 
where we have used integration by parts, (14. Il l i , ( 14. lb , Lemma |4~n Lemma|4]2]and Holder inequality. 



Jn 



Vm*) : Vii* + V(m*-^ • Vm*) : Vm*] 



r [V(m*~^ • Vm*) : Vm* + (Vm*"^ • Vm*) : Vm* + (/"^ • V)Vm* : Vm*] 
Jn 

Vm*) : Vm* + (Vm*-i • Vm*) : Vm* - div/^^Vw* : Vw* - (/^^ • V)Vm* : Vm*] 



Jn 

= "J 

Jn 



[(Vm*-i • Vm*) : Vii* + (V/-i • Vm*) : Vm* - divw^-^Vw* : Vm*] 
< C||VM*|b||V/-i|b||VM*|b, (4.23) 
where we have used integration by parts and Holder inequality. 

h = (p + A) \ [div(M*"^ • VM*)divM* + V • (m*"^ • VM*)divM*] 
Jn 

- + f [div(M*"^ • V/)divM* + divM*(VM*-^)' : Vm* + (/"^ • VdivM*)div/] 

Jn 

- (M + A) [ [div(M*"^ • V/)divM* + divM*(VM*"^)' : Vm* - divM*"MivM*divM* 

Jn 

-(m*"1 • VdivM*)divM*] 

- ifi + A) [ [(Vm*-')' : V/divii* + divM*(VM*"')' : Vm* - divM*'MivM*div/] 

Jn 
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< C||VM^||i2||V/-i||^4||V/||i4. (4.24) 
Substituting (I4.22I )- (I4.24I ) into (14.211) . and using Cauchy inequality and Corollary 14.11 we have 

^ J mW + J (/^IVm'^I^ + (ju + AMivii'^f) 

< c[OKit)f\m'-'\\i, + cm'-'wi, + cimX^ 

< C[<^K(t)f sup [(Djf(5)]2exp|c r [OK{s)]^U¥KA^)]^'ds 

0<s<T { Jo 

+C||V/-i||4, +C||V/||^,. (4.25) 

In the following, we estimate the term ||Vm'^||^4. From equation (14. Ih and (14.101 ). we know that 
satisfies 

Lw'' ^ m^ii'', in Q., 

w\x) = 0, on 5Q. ^ ' 



By ( 14.261 ) and Lemma [231 we get 

\\w%2 < C\\m^u%2 < CV^IWII V)^"'lb, 
which together with the interpolation inequality, Sobolev inequality, and ( 14. lll li yields 

< C\\VuX2\\VuXM^w\e + 

< C||Vm^||^2||Vm'^||2,[(1)^(0 + ||Vw'^||^i] 

< II V/||2, [0)^(0 + V^ll V^m'|Il2] 



< C(l)^(0l|VM'^||i2||V/||2,(l + II Vm*^M*^||i2). 
This together with Corollary 14. 11 gives 

||V/||^4 < C(l + II V^M*||i2)O^(0l|Vi<'=||2,^sup^Ojf(5)exp|c [(I)j,(^)]^[^^^^^ 

Similarly, we have 

||Vm*^-1||^4 < C{1 + \\^l^'u^-'\\:2)^K{tWu''-'\\l, ^^nv^^ 

(4.28) 

Substituting (14.271 ) and ( 14.281 ) into ( 14.251 ). we get 

m'^|M*^|2 + ^ (/^IVm^^I^ + Ou + ^)|divM^f ) 
< C<^K{t){l + \\^^u\2Wu% sup Ojf(5)exp|c r V^K{s)\'^^VVKAs)W'ds\ 

0<s<T I Jo J 

+COjf(0(l +I|V)^m*^"'|Il2)||V/-1||2, sup cD^(5)exp(c f [0^(^)]^^[4'^,.(^)] J 

0<s<T [ Jo J 

+C[OK(t)f sup [(Djf(s)]2exp|c r [OK{s)]^Ui'K,M]^ds\. 

0<s<T { Jo ) 

Integrating over (0, t) and using ( 11.61 ). we complete the proof of Lemma 14.41 
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Note that T > and r e (3, 4] are arbitrary for all the above estimates which will be useful to get 
the blow-up criterion of the solution in the next section. To obtain the strong solutions, we have to 
take T small enough. Therefore, we assume T e (0, 1). Moreover, we take r = 4 for simplicity. 

Suppose foT I < k < K 



Wu'" ^\\lHo,T;W^-^) + \\u'' 'llL~(0,r;//2)+ sup f m'' 'P<Mi, 

0<t<T JSi 



(4.29) 



for Ml > 1 large enough. 



Throughout the rest of the section, we denote by C a generic positive constant which may be 
dependent of fi, A, Q., mo, «o> "o and other known constants but independent of M\, K, 5 and T. 



Lemma 4.5. Under the conditions ofTheorem \4.1\ we have for any k > I and T < M^^ 

' sup J^mVfP + Jo^ i IVwfP < CM\, 

Q<t<T 

\\m%«'(Q,T;W''i) + ll"*'llL»(o,r;W"-'?) ^ C, 
^ < ^Qm*" <n'' < lorn'', in Qj, 

Proof. By ( 14.291) . (I3l4l) . (1411) i -(1411)9 and Sobolev inequality, we have 



' sup ||m*(Ollwi,, < Cexp{CMiT2 

0<t<T 

sup \\n\t)\\wi, < Cexp{CMiT-2] 



\Q<t<T 



Denote Ti = M'^, we get for T < Ti 



sup llm'^COIIv^i,. < C, 

0<f<7" 

sup \\n''{t)\yu, < C. 
iO<t<T 



(I4.30l) i and Sobolev inequality give 



This together with (14.51 ) implies 



sup O/f(0 < C. 

o<f<r 



By dlH), (US), (14301) and (1432]) . we have 

We obtain from ( 14.291 ). (14.311 ) and Sobolev inequality 

sup ^^,4(0 < CMf . 
o<f<r 

It follows from Corollary KH <S2B, and K34\ that for T < 

IIV/|lL2(o,r;L6) + IIV/|lL»(o,r;L2) < Cexp {CTM 



(4.30) 

(4.31) 
(4.32) 

(4.33) 
(4.34) 
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Take T2 = M^^ we have for T < T2 



I|Vm*IIl2(o,T;L«) + IIV«^llL»(0,r;L2) ^ C. 



(4.35) 



By Lemma|431 (l429l) . (g^B, (l434b . (l435]) . and Sobolev inequality, we get for < f < T < T2 

< C + £ Igp + C j^' (1 + II V^^M^^II^z) ||Vm*^||2, + c (1 + II V)^m'"'|Il2) ||Vm^-i 11^ 

< C+ r |gp + c(l+ sup II Vm^M^ 1 1^2 1 + CMj^r. 



Using Cauchy inequality, we have for T < T2 



sup r m*iM*p+ r r ivm*^p<c+c r i^i^ +cMfr. 

o<f<r Jn Jo Jn Jn 

By (I4J91) . LemmaO (l43n) . (14331) and (l436l) . we get for T < 

+ ll«'llL~(0,r;//2) ^ C + C r Igp + CMfr. 

Jn 

Since T < M~*^, we have 

sup r m'=iM'=|2+ r r ivM*=|2+iiM'=ii^2(or.,^2,,)+iiM^ii^„(o^.^2)<c+c r i^p. 

o<r<r Jn Jo Jn Jn 

Let Ml > C + C |gp, we obtain for T <T2 

sup r m'=|M'=|2+ r r |VM*^|2 + ||/||^2(o,r;W2.';) + ll"'^llL»(0,r;//2)<^l- 

^Jn Jo Jn 



(4.36) 



(4.37) 



0<t<T ■ 



By induction, (14.371 ) is valid for any k € [1, ^]. Since Mi is independent of K, and K is arbitrary, we 
conclude that (14.371 ) is actually valid for all ^ > 1. 
From (14.371 ). we obtain for T < T2 



sup r mv?p+ r r ivMf|2< cm^. 

o<r<r Jn Jo Jn 

Summarily, we have for any ^ > 1 and T < T2 

'sup X,mW + /;/jV.fp<CMf, 

0<t<T 

\\'^%HO,T-W^'i) + ll"''llL"(0,r;ffln//2) ^ ^1' 

IN*'llL~(0,r;»'i.'?) + Wn'^WL'-iOJ-W'") ^ 
< Sj^ni' <n^ < Tom*', in Qj, 

{\\P'\\L-iQr) + \\P'j\L-iQr) + iKh^iOr) < C 

The proof of Lemma 4.5 is completed. 



(4.38) 
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We will show that the full sequence (m'^, n!^, u^) converges to a solution of (ll.ll )- (ll.5l ). To do this, 
we denote 

7J7^+i = ni'^^^ _ ,n\ n'^' = n^^' - n\ u'^' = u'^' - and p'^' - P'^' - pK 
It follows from (14. Ih that 

m'^+iMf+i + m*^+i/ • V/+^ - Lm*^+' + VP*'^^ = m''^\-u'^ - u'' ■ V/) - m'^u'' ■ VuK (4.39) 
Multiplying (14.391 ) by u'^'*'^, and using (14.381 ) and Sobolev inequality, we have 

2 <i? Jn Jn 

= r p'^^'divIZ^^+U r m'^+i(-Mf-/.VM'=)77*^+i - r mV-VM*^)-n'=+' 
Jn Jn 

< ||p'"''||i2||diVM*^+l||^2 + \\m'-'\2\\u'} + i/*^ • V/||^3||i7'=+l|b + II ^'U%2\\ V)^||L~l|n'^^'lbl|V/||i3 

< ||P'^' 11^2 ||divM*^+' 11^2 + C||m'^+'||i2(||VMf||i2 + M^,mu'*'h2 + C\\ ^^^^u\2\m'^\2M, 

^ Jn 
This gives 

- r m^+i|I7'^+V+ r [;/|Vn'^+f + (// + ^)|divI7*^+ip] 
Jn Jn 

< C||n'+i||2, + C||m'^i||2,(||VMf||2, + M\) + CM^II 4^^%. (4.40) 
From (I4.1l )i . we have 

mf*^ + m^^-^MivM^^ + m^divM^ + ■ 11^^^ + ■ Vm^ - 0. (4.41) 
Multiplying ( 14.411 ) by nt''^ , integrating over n, and using integration by parts and (14.381 ). we have 

1 r i^'^+ip = _ r |5H'=+i|2divM'^-2 r (m^^divM^^ + M^^- Vm*^)m^+i 
Jn Jn Jn 

< C||/||^2,,||m'^+i||2, + C||m'=+'||i2(||VI7^||i2 + ||I7*^ • Vm^^ll^j) 

< C||M*^||^2,,||m'^+i||2 + C||m'=+'||i2(||Vi7*^||i2 + ||Vm*^||^3||M'^||i6) 



Hi: 

Similarly, we have from (I4.1l i7 



< C||/||,^2,,||m'^+'||2 +C||m'=+'||i2||VM*^||^2. (4.42) 



^ £ IH'^^f < C||/||iy2,,||H*^+^||2, + C||H*^+>|L2||Vm'^||^2. (4.43) 

By (I4.42I) - (I4.43I) and Cauchy inequality, we have 

^ J^(l^'^¥ + l«'^¥) < C||/||,^2,,(||m*^+i||2, + IIH^^+^II^,) + C(||m^-^^||i2 + ||H'=+1||^2)||Vm*^||^2. (4.44) 

By (14.401) . (14.441 ) and Cauchy inequahty, we get 

^ r (|m^-+'p + |H'=+'|2 + m'^+i|I7'=-^ih+ f MVI7^-+'|2 + (;, + ^)|divI7*^-^ip] 
Jn Jn 
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£ Jn 

for any e > 0. Denote 

Ui\Wu^^'\^ + (ju + A)\divu''^'fm, 





- k 






^'^\t) -- 








D\s,t) -- 





11^2 ■ 



We have 

This together with Gronwall inequality, (I4.38l l and ^^'^"'^'(0) - implies 

ip''*\t) + >^''^\s) < exp (c + ^) [£i//''(s) + CM2/(5)], (4.45) 

where C depends on M\ and other known constants related to C. 
By (14.451) . we have 



sup /+H0 + r (/'''"'k'^) < exp (c + — ) (CM^r + e) ( sup /(t) + f tfr^s) 

0<t<T Jo \ £ / \0<t<T Jo 

Take e = — exp (c + c) {CM\ + 1), and To = min{r2, e), we have for T = Tq 
16 ^ ' 

sup [/+ko]^ +( r V'^'(^)]' < i sup (At)f + i f 

<t<To \Jo I ^ 0<t<To \ Jo 



0<f<7'o 

This implies 



y j sup [<p'^'{t)]-2 i/f'-^'is)]' \ < cx,. 

ti \ o<t<To \Jo / 



k=l 



Recalling the notations of <p (t) and ijj {t), we get 

00 

(ll^'"'llz.»(0,ro;L^) + IIH'"'llL»(o,ro;L2) + ll^'"'llL^(0,ro;//i)) < ~- (4-46) 



/t=l 



Denote m - ^ + m\ « = ^ + n^, m = ^ m', we have 

(=2 i=2 ;=1 



IN* - m||i,»(o,ro;L2) + II"*' - «llL"(o,ro;i^) + H"^ " "Hl^cojo;//^ 
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CO 




CO 




CO 




+ 




+ 




i=k+l 




i=k+l 


L»(0,ro;i2) 


i=k+l 



(4.47) 



< > I IN'llL-co.ro;!^) + ll«'llL»(o.ro;L2) + ll"'llL2(o,ro;ffi) ) ^ 0' 
i=k+i 

as k ^ oo. Here we have used (I4.46I ). 

Therefore, we conclude the convergence of the full sequence (m*, u'^) as k oo 

'm^^m, in L°°(0, Tq; L^), 
n^^n, in L°°{0,Tq;L^), 
u'^^u, in L^{Q,Tq;HI). 

It follows from Lemma 1431 (14.471 ) and lower semi-continuity of norms, we conclude that (m, n, u) is a 
solution to ( |l.l| )-( fT31 ) under the conditions of Theorem 14. 1[ and that the following estimates uniform 
for 6 are obtained 

{ sup LmkP + /;°/^|V«,p<C, 

0<t<TQ 

M\lHo,To-w^'i) + ll"llL~(o,ro;//>ff2) + ll'«||i,oo(o^7.g.iyi.,) + WnWu^^Qj^.y^u,) < C, 
^ < < « < Tom, in Qt^, 

\\P\\l-^(Qto) + ll^mllL"(ero) + ll^nllL^lCro) ^ 

The uniqueness can be obtained similar to the proceeding of the convergence of full sequence. The 
proof of Theorem 14. 1 l is completed. 



Proof of Theorem ll.il 

1q = niQ + 6, ?1q 



Denote mf, = niQ + 6, nf. = jiq + 6, we have as 5 ^ 



(m^ ^ mo, in W^'^, 
\nf.^no, in W^-^. 



(4.48) 



Since < ^gmo < no ^ soMq, in Q, without loss of generality, we assume < < 1 < sq, we have 



By Lemma [231 we can find a e D for each 5 > such that 



(4.49) 



Lul + VP(m^, n^) = ^m^g, in Q, 
I "ol'5n ^ 0. 



(4.50) 



It follows from ( 1448] ). (|430l ). and Lemma|23]that 
as 5 ^ 0. 

Consider (ll.ll )- (ll.5l ) with initial data (mo, no, "o) replaced by (m^, n^, m^), we obtain from Theo- 
rem |4?T]that there exists a To > independent of 6 and a unique solution {m^ , n^, u^) for each 5 > 
with the following estimates: 

f sup X,'^W + /o'7nl^"fl'^C' 

o</<ro 

Wu^HOJo-W^-i) + ll"^llL»(0,ro;//>ff2) + W^^^imuW^") + ll"''llL"(0,ro;H"-'/) ^ C, 

^ < £o'«'' ^ ^ ^o'w'', in Qto, 
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where C is a positive constant, independent of d. 

Taking 6^0 (take subsequence if necessary), and using the similar arguments in IH, under 
the conditions of Theorem 11.11 we get a solution (m, n, u) to (ll.ll )- (ll.5l) with the regularities like in 
Theorem ll.il The uniqueness can be proved by the similar arguments in the proof of Theorem 14. 1 1 
The proof of Theorem ll.ll is completed. 



Let {m, n, u) be a strong solution to the problem (ll.ll i- (ll.5l ) in Qj with the regularity stated in Theorem 
11.11 We assume that the opposite holds, i.e. 



In this section, we denote by C a generic positive constant which may depend on A, Q, mo, hq, uq, 
M, T* , and the parameters in the expression of P in (11.51 ). 
Similar to Lemma l4~n we get the first lemma: 

Lemma 5.1. Under the conditions of Theorem U .2\ we have for all < T < T* 



o<t<T Jn 

Proof. Multiplying il.lh by r\u\''~^u, integrating the resulting equation over O, and using integration 
by parts, we obtain 



5 Proof of Theorem fL2 



Um sup ||m||L~(o,r;L~) < M < oo. 



(5.1) 



Sj^m <n< SQin, in Qj. 
Lemma 5.2. Under the conditions of Theorem U .2\ there exists some r € (3,4] such that 





(5.2) 



For ei € (0, 1), define 




Case 1: If 




(5.3) 



A direct calculation gives 




(5.4) 



and 




(5.5) 



By (|5^ and (1531) . we get 
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= rj div(\ur^u)P-r{r-2)(ii + A) J ImI'-^m • V|M|div 



r(r - 2)(ji + A) 



f \ur 

Jn 



^\u-V\u\f 



(5.6) 



where we have used Cauchy inequality. By (15.41 ) and (15.61 ). we get 



< r 



This together with ( 15.11 ). ( 15.31 ) and Cauchy inequality implies for e > 



< r 



< C 



r ii/r~'|viM 

Jn 



+ <f>i£i,r)\ fir 



[ \ur^\^\u\\\ 

Jn 



where we have used P{m, n) < Cm2 , which can be obtained from (15.11 ). Lemma l5?T] and the expression 
of P. 

This together with ( 15.31 ) and ( 15.41 ) gives 



/i(r- l)-0(ei,r) // 



Jn 

£[1 + m,r)] £ iMr-2|viMi|' + ^ 



n 

Taking £ ^ [1 + 0(£i, r)]-ir - 1) - 0(gi , r) ( ^^'"^ -//)], we get 



* Jn 



C[l+0(£i,r)] 



4//r(r - 1) - 0(£i , r) [3r(r - 2)(// + /I) - 4//r] 



/ 

Jn 



m\u\ 



(5.7) 



for any r € (3,4]. 
Case 2: 



By (l52l ). we get 

^ J" m|Mr + J" r|Mr-2|^|yj^|2^(-^+^)|div„|2+^(;,_2)|V|Mf) 



(5.8) 



^ r I div(|Mr"^M)P - r(r - 2)(jU + /t) I divM|M|"^|M|"^M • V|m| 
Jn Jn 



Ar-2, 



Jn 



/■-2 I--4 



20 



< c 



where we have used Cauchy inequality, (I5.1l l and P < Cm5. Therefore, 

— f m|Mr+ r /irlwr^lVMp + Mr - 2)r f |Mr-2|v|M|p 

< c r,„^i.nv.i.^l^:^^^^ r N-|vN|l (5.9) 

Jn ^ Jn 

It follows from (I5.4I ). (15.91 ) and Cauchy inequality that for any eo € (0, 1) 

— 

. c J^.s?M-|vM|.,.../j„r|v(^)r.^(J^„,i„r) ' 

^ Jn 
Combining (15.81) . we have 

d r , r (r-2)2(u + i)l f ,1 |2 
- J^'^l«r + '- Ml-£o)0(ei,r)+M'--l)-^ -f- J^M'''|VM| 

I--2 

< C r m'sr |m|'-2| v|m|| + f ml^r ) ' . (5.10) 
Jn ' ' 4//reo \Jn / 

(Sub-Case 2i): If 3 € (2 + 00), i.e. n < 3/1, we have for r € [3, 00) 

AsiLiir - 1) 

cf>{£ur) = — , (5.11) 

^ 3ir - 2)(jj + A) - 4iu ^ ' 

Define ^ 

/(£(), ei,r) ^ ;u(l - E(i)<p{£\,f) + ^{r - 1) - ^^^^ — 2)^(// + /I) ^ ^^^^^ 

By (ISTTTl) and (I5l2l ). we have 

(r-2)20u + .t) 

/(eo,ei,'') = ^iiir-\) , (5.13) 

3(r - 2)(/i + /I) - 4yu 4 



for r e [3, 00). Particularly, 



/(o.i,3)=Ji;l.Zii^>o. 

3/1 4 



Here we have used ^ < X <'^[i. 

Since /(eo. £i> ?") is continuous w.r.t. (sq, s\,r) in (0, 1) x (0, 1) x [3, 00), there exists eo. £1 £ (0, 1) 
and r\ € (3, 4] such that 

/(eo,ei,ri) > 0. 
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From (IS.lOl l. Cauchy inequality and Holder inequality, we obtain 
j^J m\up +nf(so,Eun) J \ur~^\W\u\\^ 

rj-2 rj-2 

< nf{£o,£i,n) f \ur-^\V\uf + ——^ -If mluA" +— mluA" . 

Ja 4n/(eo,ei,n) \Jn / 4//ri£o \Jn / 

Therefore, 

^ f m\ur<C[—-^ -, + —^if MuA^' . (5.14) 

(Sub-Case 22): If 3 ^ (2 + cx,), i.e. fi > 3A. 

In this case, it is easy to verify that the following inequality holds for any re (3,4] 



//(I -Eo)(f>{£ur)+fi{r- 1) 

We obtain from (15.101 ). (15.151 ). Cauchy inequality and Holder inequality 

|2 



> 2iu. (5.15) 



< C r m'5F|M|'-2|V|M|| + — ^1 r m\u\' 



1--2 



r-2 



This implies 



- r m|Mr<— I r m|Mr| ' . (5.16) 
dt Jn IJn / 



Particularly, (15.161 ) is also valid for r = ri . 

Summarily, for Case 2, we obtain from (15.141 ) and (15.161 ) 



m|Mr'<c|j" m|M|'''J '' , (5.17) 



for some constant C,if A < and (15.81) is valid. 
It follows from (15^ and (15.171 ) that 



n-2 

m|Mr'<c|j" m|M|'''J , (5.18) 



for some constant C,if A < ^fi. 



3 ' 

Since e (0, 1), we complete the proof of Lemma |5]2] after using young inequality and Gron- 
wall inequality in (15.181) and still denoting ri by r. □ 

Similar to Lemma l4!2l and Corollary 14.11 we get the next lemma. 
Lemma 5.3. Under the conditions of Theorem U .2\ and ( 15.71 ). for < T < T*, we have 

i\\P\\L-iQr) + \\Pm\\L-iQr) + WPnh-iQr) < C, 
\l|VM|lL2((),r;L6) + l|Vi<llL»(0,r;L2) ^ C. 
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Here we have used (15.11 ) and Lemma [5^ 
Lemma 5.4. Under the conditions of Theorem U .2\ and ( 15. il ), for < T < T*, we have 



r miiip + r r iv^p < c, 

Jn Jo Jn 



where ii - Uf + u ■ Vu. 

Proof. Similar to Lemma l4!4l we have 



r m|Mp+ r r (/i|VM|^ + Ou + i)|divMp) 
Jn Jo Jfi 

' j| [(l + ||V^M|li2)||VM||2,] 

< C + cj^ (II V?^M||2, II Vm||2,), 



< C + C 



where we have used (15.191 1 and Cauchy inequality. This together with Gronwall inequality and (15.19) 1 
completes the proof of Lemma 5.4. n 



Denote w = u - h, where h is the solution to 

(Lh = VP{m, n), in Q X (0, T], 
\h\an = 0. 

Similar to (I4.26I ). we have 

Lw = mil, in Q.x{0,T], 



(5.20) 



w\dii = 0. 

Due to (15.11) . (15.211 ). Lemma [231 and Lemma |5^ we immediately give the following result. 
Corollary 5.1. Under the conditions of Theorem \1.2\ and ( 15.71 ), for < T < T*, we have 

\M\l^(0,T;W^.i) ^ C. 

In the following, we give the estimates of the derivatives of m and n. 
Lemma 5.5. Under the conditions of Theorem U .2\ and ( 15. iP . for < T < T*, we have 

sup ||(Vm, V?i)(OIIl'' < C. 

te[0,T] 

Proof. Differentiating the equation (ll.ll )i with respect to Xi, then multiplying both sides of the result- 
ing equation by ^|5,m|^"^5;m, we get 

dt\dim\'' + div(|<9,m|''M) + {q - l)|<9,m|'^divM 

+qm\dim\'^~^dimdidivu + q\dim\^~^dimdiU ■ Vm = 0. (5.22) 
Integrating (15.221 ) over Q, we obtain 

— r |Vm|'? < C \ IVwIIVml'' + <7 f m|VdivM||Vm|«-i 
dt Jn JQ. Jn 

< C||VM||i»||Vm|||,+C||V2M||i,||Vm|||;^ (5.23) 
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Similarly, 

I £ \vnf < c\m\L'^\m\i, + cw^/UMmf-'. (5.24) 

By S5?23\i - K24\i . we have 

;^ (llVmlli, + llVnlli,) 
at 

< C(l + ||Vw||l» + \\Vhh^){\\Vm\\L, + ||V«||i,) + QIVMIl" + QlV^hh.. (5.25) 
Since P,„ and P„ are bounded, we have from (15.201 ) and Lemma 1231 

||V2/j||i, < C(||Vm||i, + llVnlli.). (5.26) 
Applying (I5.25I )- (I5.26I ) and Sobolev inequaUty, we get 

^ (llVmlli, + \\Vn\y) 
at 

< C(l + ||w||iy2., + ||V/j||L»)(||Vm||L, + ||V«||l,) + C\\vMl^. (5.27) 
Using (I5TT91) . (l5^ . Lemmas |23]|231 we have 

\m\L- < c(l + \m\BMOin)Me + \\^^h\\L.)) 

< C(l + ||P|L»nL2 ln(e + ||VP||iO) 

< C(l + ln(e + llVmlb + ||V«|b)). (5.28) 
From (I5.271 l- (l5.281 l and Cauchy inequality, we get 

4 (IIVmllL, + ||V«|k,) 

< C(l + ||w||^2,, + ln{e + \\Vm\\L, + ||Vn||„)) (||Vm|k, + WVnh,) + C\\VMu'- (5.29) 
Denote G{t) = e + ||Vm(OllL« + I|V«(OIIl'?, we have from (15.291 ) 

4g(0 < CIIV^wIIl, + C(l + \\w\\w2.,)G{t) + CG(0 lnG(0. (5.30) 
at 

Multiplying (15.301 ) by and using G > 1, we have 

^lnG(0 < C(l + ||w||iy2,,) + ClnG(0. (5.31) 

Using Gronwall inequality, Corollary 15. H and (15.31| l. we complete the proof of Lemma 5.5. □ 
Lemma 5.6. Under the conditions ofTheorem \1.2\ and ( 15. i I ), for < T < T*, we have 

Proof. Rewrite (11.11 )^ as 

Lu - mil + VP{m,n). 
By (15. lb . ( 15. 19b . Lemma 1231 Lemma [53] and Lemma [531 we have 

\\u\\h2 < C(|M||^2 + llVPIb) 
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< C{\\m^u\\i2+\\^m\\i2+\\^n\y)<C, 



and 



M\l\q,t-w^.i) ^ C(||mM||i2(o,r;L'?) + \\^ P\\lHo,t-ui)) 

< C(||VM||i,2(oj.i2) + ||Vm||i2(o,T;L'?) + I|V«|Il2(0,T;L9)) ^ C. 



□ 



By ( I5.1I ). Lemma Lemma 1531 and Sobolev inequality, we get the following result. 



Corollary 5.2. Under the conditions of Theorem \1.2\ and d5.iD . for < T < T*, we have 



n Jo Jn 



m\utf+ \Vutf<C. 




By (I5.1I ). Lemma I5.3[ Lemma 15.51 Lemma 15.61 and Corollary 15.21 we know that T* is not the 
maximal existence time for the strong solution (m, n, u){x, t) to the problem (|1.1| )- (|1.5I ). This is a 
contradiction with the definition of T* . Therefore, (15.11) is invalid, i.e. 
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Abstract 

In this work, we prove the existence and uniqueness of the global classical spherically 
symmetric solution to the compressible isentropic Navier-Stokes equations with vacuum 
in a bounded domain or exterior domain il of M"(n > 2). This is a extensive work of 0, 
where the global classical solution of compressible Navier-Stokes equations in one dimen- 
sion was obtained. As pointed out in |6j, the regularity of velocity u e Hl^^{[0,oo); H^) 
could not be improved to C^{[0,oo);H^), otherwise, it would blow up in finite time. 
Compared to [B] , the regularities of velocity with respect to the space variables are im- 
proved. The analysis is based on some new mathematical techniques and some new useful 
estimates. This can be viewed to be the first result on global classical solution with large 
initial data and vacuum. 

Keywords: Compressible Navier-Stokes equations, vacuum, global classical solution. 
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1 Introduction 



In this paper, we consider the initial-boundary value problem of compressible isentropic 
Navier-Stokes equations in a bounded domain or exterior domain 0, of M"(n > 2): 



(1.1) 



Pt + V- (pu) = 0, p > 0, 

{pu)t + V • (pu ® u) + VP(p) = ^Au +{fi + A)V(V • u) + pi, 

for (x, t) £Qx (0,-|-oo), where 

O = {x a < |x| < 6}, < a < 6 < oo, f(x,t) = /(|x|,t) — , 

|x| 

p and u = (ni,n2,--- ,n„) denote the density and the velocity respectively; P{p) = Kp^, 
for some constants 7 > 1 and K > 0, is the pressure function; f is the external force; the 
viscosity coefficients p and A satisfy the natural physical restrictions: p > and 2p + nX > 0. 
We consider the initial condition: 

(P, ^)\t=o = (Po, uo) in n, (1.2) 

and the boundary condition: 

u(x, t) — )• 0, as |x| — )• a or 6, for t>0, (1-3) 

where 

Po(x) = Po(|x|), uo(x) = no(|x|) — . 

|x| 

We are looking for a classical spherically symmetric solution (p, u): 

p{x,t) = p{r,t), u{x,t) = u{r,t)-, 

r 

where r = |x|, and {p,u){r,t) satisfies 

pu 

pt + [pu)r + m— = 0, p > 0, 

[pu)t + (pn )r + m—- + Pr = l^[Ur + "I-),. + p/, 

for (r, t) G {a,b) x (0,oo), with the initial condition: 

{p{r,t), u{r,t))\^^Q = {po{r), Uo{r)) in I, (1.5) 

and the boundary condition: 

n(r, t) — )• 0, as r — )• a or 6, for t > 0, (1-6) 

where m = n — 1, i/ = 2p-|-A> p > and / = [a, b]. 

Let's review some previous work in this direction. When the viscosity coefficient p is con- 
stant, the local classical solution of non-isentropic Navier-Stokes equations in Holder spaces 
was obtained by Tani in [2D] with po bounded below away from zero. Using delicate energy 
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methods in Sobolev spaces, Matsumura and Nishida showed in |17[ [T8] that the existence of 
the global classical solution provided that the initial data was small in some sense and away 
from vacuum. There are also some results about the existence of global strong solution to the 
Navier-Stokes equations with constant viscosity coefficient when po > 0, refer for instance 
to [U [12] for the isentropic flow. Jiang in [12] proved the global existence of spherically 
symmetric smooth solutions in Holder spaces to the equations of a viscous polytropic ideal 
gas in the domain exterior to a ball in M" (n = 2 or 3) when po > 0. For general initial 
data, Kawohl in [T^ got the global classical solution with po > and the viscosity coefficient 
H = fi{p) satisfying 

< Po < p(p) < Pi, for p > 0, 

where po and pi are constants. Indeed, such a condition includes the case p(p) =const. 

In the presence of vacuum. Lions in [16] used the weak convergence method to show 
the existence of global weak solution to the Navier-Stokes equations for isentropic flow with 
general initial data and 7 > | in three dimensional space. Later, the restriction on 7 was 
relaxed by Feireisl, et al [9] to 7 > |. Unfortunately, this assumption excludes for example 
the interesting case 7 = 1.4 (air, et al). Jiang and Zhang relaxed the condition to 7 > 1 in [13] 
when they considered the global spherically symmetric weak solution. It worths mentioning 
a result due to Hoff in [llj, the existence of a weak solution for positive initial density has 
been proved when 7 = 1. 

There were few results about strong solution when the initial density may vanish until 
Salvi and Straskraba in [19J, where 17 is a bounded domain, P = P{-) £ C^[0,oo), po G H^, 
uq S Hq f]H^, and satisfied the compatibility condition: 

Luo(x) - VP(po)(x) = VPog> for g e L\ (1.7) 

Afterwards, Cho, Choe and Kim in [2l[3lll] established some local and global existence results 
about strong solution in bounded or unbounded domain with initial data different from [19] 
still satisfying (|1.7p . Particularly, Choe and Kim in ^ showed the radially symmetric strong 
solution existed globally in time for 7 > 2 in annular domain. As it is pointed out in [1] 
that the results have been proved only for annular domain and cannot be extended to a ball 
Q = Br = {x G : |x| < i? < 00}, because of a counter-example of Weigant |21j . Precisely, 
for 1 < 7 < 1 -|- Weigant constructed a radially symmetric strong solution (p, u) in 

Br X [0,1) such that \\p{-,t)\\L°°{Bji) — 00 as t — )■ 1^. A recent paper ^ written by Fan, 
Jiang and Ni improved the result in [4J to the case 7 > 1. 

The local classical solution was obtained by Cho and Kim in [5] when the initial density 
may vanish and satisfying the following compatible conditions: 

Luo(x) - VP(po)(x) = po[gi(x) - f(x,0)], (1.8) 

for x G 0, gi G Dq and ^/pogl G L^- Recently, we used some new estimates to get a unique 
globally classical solution p G C^{[0,(yo); H^) and u G HI^^{[0,og); H^) to one dimensional 
compressible Navier-Stokes equations in a bounded domain when the initial density may 
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vanish, cf. [6j. Since Xin in [22j showed that the smooth solution {p, u) G C"'^([0, oo); /f'^(M^)) 
must blow up when the initial density is of nontrivial compact support, so the regularities of 
u with respect to time variable obtained in [6j could not be improved, refer to [6j for more 
details. 

Since the system (|1.4|) have the one dimensional feature, the results in [B] are possible to 
obtain here. Besides, the regularities of u with respect to space variable could be improved 
though these of u with respect to time variable couldn't be done. This causes some new 
challenges, which are handled by some new estimates. 

This can be viewed to be the first result on global classical solution with large initial data 
and vacuum. 

Notations: 

(1) Qt = Ix [0,T], Qt = TIx [0,T] for r>0. 

(2) For p > 1, LP = LP{0,) denotes the space with the norm || • \\lp- For k > 1 and 
p > 1, W'^'P = W^'P{Q) denotes the Sobolev space, whose norm is denoted as || • H^yfe^p; 

ffk ^ p^fe,2(j^)_ 

(3) For an integer k > and < a < 1, let C'^^"(r2) denote the Schauder space of 
function on fi, whose A;th order derivative is Holder continuous with exponent a, with the 
norm || • ||(^fc+ti. 

(4) For an integer A; > 0, denote 



(7) L := fj.A + {fi + A)Vdiv is the Lame operator. 
Our main results are stated as follows. 

Theorem 1.1 Assume that po > satisfies po £ H H'^,Pq € H'^, uq G n Dq and 
f G C{[0,oo); H^), fj G L^^^([0, oo); L^), and the initial data pQ, Uq satisfy the compatible 
condition il.8\) with gi(x) = gi{r)j. Then for any T > 0, there exists a unique global 
classical solution (p, u) to il.l\) - [T73\) satisfying 




with the norm 




(5) D'''P={v G Ll^{n)\\\V''v\\LP <oo},D'' = D'''^. 

(6) = Dp'^ is the closure of C^{n) in D^^^. 



{p,p^)eCi[0,T];H^), p>0, ipu)tGCi[0,T];H'), 
uGC([0,r];Z?3nZ)i), ut e L^{[0,T];Dl)nL^{[0,T];D^), e L^{Qt). 
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Remark 1.1 (i) Note that if is bounded and locally Lipschitzian, then D^'^ = W^'^. See 
|10j for the proof. 

(ii) By Sobolev embedding theorems, we have 

H''{I) ^ C^'^I), for A; = 1,2,3, 

this together with the regularities of (p, u) give 

{p,p^) G C([0,r];Ci+l(/)), u G C([0,r];C2+l(/)). 
Since (p(x, t), u(x, t)) = {p{r,t),u{r,t)f), we get 

ip,p'') G Ci[0,T];C'+Hn)), u G C7([0,T];C2+5(n)), 
which means (/O, u) is the classical solution to (ll.ip - (ll.3p . 

Theorem 1.2 Consider the same assumptions as in Theorem \1. 11 and in addition assume 
that po G H^pI G i^^ V(^) G ^i^V^gi G L^, poV^gi G L^, uq G D^, f G 

C{[0,ooy,H^)nLl^{[0,^y,H^), it G C([0,oo);/7i)nL2^,([0,oo);/72) ^^^f^^ ^ L2^^([0, oo); L^). 
T/ien the regularities of the solution obtained in Theorem can be improved as follows: 

ip,p^)GCi[0,T];H^), ^GW^'^iQr), ipu)t G L^{[0,T]; H^) n L\[0,T]- H''), 
uGL~([0,r];i?ini?5)nL2([0,T];D6)^ G L~([0, T]; Z)i) n L2([0, T]; I)^)^ 

(VpV2ut,pV3ut) G L-([0,r];L2), ^V^u^ G ^^(Qr), P^^tt G ^'(Qt), /O^u^t^ G L\Qt), 
pluuGL°-{[0,T];L^)nL\[0,T];H^), p^uu G L°-{[0,T]; H^) D L\[0,T]; D^). 

The rest of the paper is organized as follows. In Section 2, we prove Theorem 11.11 In 
Section 3, we prove Theorem 11.21 by giving some estimates similar to and some new 
estimates, such as Lemma 13.61 Lemma 13.71 and Lemma 13.81 

The constants u and K play no role in the analysis, so we assume ly = K = 1 without 
loss of generality. 

2 Proof of Theorem 11.11 

In this section, we get a unique global classical solution to (jl.4p - (jl.6p with initial density 
Po ^ ^ > and & < oo by some a priori estimates globally in time based on the local 
solution. Moreover, the estimates are independent of b and 6. Next, we construct a sequence 
of approximate solutions to (ll.4jl - ()1.6p under the assumption po ^ > 0. We obtain the 
global classical solution to (|1.4p - (jl.6p for po > and b < oo after taking the limits (5 — )• 0. 
Based on the global classical solution for the case of 6 < oo, where the estimates are uniform 
for b, then we can get the solution in the exterior domain by using the similar arguments as 
that in [3]. 

In the section, we denote by "c" the generic constant depending on a, 1 1 po 1 1 ) 1 1 Po 1 1 ) 1 1 ^o 1 1 , 
||uo||£)3, T and some other known constants but independent of 6 and b. 
Before proving Theorem I l.H we give the following auxiliary theorem. 
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Theorem 2.1 Consider the same assumptions as in Theorem and in addition assume 
that po > 6 > and 6 < oo. Then for any T > 0, there exists a unique global classical 
solution {p,u) to l^l.4\ )- [T^) satisfying 

peC{[0,T];H\l)), p>|, n«GL2([0,T];L2(/)), 
ueC{[0,T];HHl)nHUl)), ut€C{[0,T];H',{I))nL\[0,T];H\l)). 

The local solution in Theorem 12. II can be obtained by the successive approximations like in 
O [5], we omit it here for simplicity. The regularities guarantee the uniqueness (refer for 
instance to [21 ^3]). Based on it, Theorem 12.11 can be proved by some a priori estimates 
globally in time. 

For T G (0, +oo), let {p, u) be the classical solution of (ll.4p - (|1.6p as in Theorem 12. II Then 
we have the following estimates (cf. [4] and [8]): 

Lemma 2.1 For any < t < T , it holds 

\\{p,p')\\H? + \\{Pt,{p')t)\\m^+ r\\{pu,{p')u)\\h<c, p>-, 

Jo c 

and 

[ {r'^pul + r'^ul^ + r'^ul + r™-\2) + f {r^ul^ + r'^-^u'i + r™^/^^,) < c, 
J I Jqt 

where /i2)||x = ll^illx + ||/i2||x5 for some Banach space X, and hi £ X, i=l,2. 

Remark 2.1 {i) For the estimates about p^ , since p and p"' satisfy the linear transport 
equations: pt + u- Vp + pV • u = and {p'^)t + u • V{p^) + jp'^V • u = respectively, then by 
using the similar arguments as that of Lemma 3. 6 in [4^, we get 

III \V^pix,t)\^dx + I |v2(pT)(x,t)|2dx| 

< c||Vu(.,t)||^. (||Vp(-,t)||i. + ||V(/,^)(-,t)||^0 
+c{\\V'G{;t)\\L^ + l|V'(p^)(-,t)||L.) {\\V'pi;t)\\L^ + \\V\p'r){.,t)\\L^) 

< c{\\Vp{;t)\\m + \Nip'){-Mm + \\G{-Mh2) , 

where G = i^V -u — p"' is the effective viscous flux; we have used the estimates ||u(-,t)||£,i < c, 
\\u{-,t)\\£,2 < c and \\p{-,t)\\fji < c in Since \\G{-,t)\\f^2 £ L^(0,T) given in 14], hence, 
an application of the Gronwall inequality gives 

\\{p,p'r){.,t)\\H2<c, (*) 

where we have used the following Sobolev inequalities for radially symmetric functions defined 
in ri' = {x G M"; |x| > a > 0}.- ||V/9||loo < c|| V/9||j:^i , ||Vu||loo < c||Vu||j:^i. From (*) and a 
direct calculation, we get 

\\{p,P^)\\m^<c. 

[a) For the case a > and 7 different from that in J 21], it is interesting to investigate 
the existence of global strong or classical spherically symmetric solution. We leave it as a 
forthcoming paper. 
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From (jl.4p . we get 

Ur H + pf. 

r / r 

Differentiating (|2.ip with respect to t, we have 



(2.1) 



pUtt + PtUt + PtUUr + PUfUr + PUUrt + ip'^)rt = Urrt + mr '^{rUrt - Ut) + Ptf + pft- (2.2) 

Lemma 2.2 For any < t < T , it holds 

J I JQt 

Proof. Multiplying (j2.2p by r^uu, integrating over /, and using integration by parts, Lemma 
12.11 and Cauchy inequality, we have 

, ""/^^^ + \jt /(""^'t + "^^'""'^?) 

r"^{Ptf + pft)utt - / r'^iPtUt + PtUUr + putUr + pnn,.t + {p'^)rt]utt 
I Jl 

< - r'^ptututt - r'^iptuur - ptf)utt + c r"^ pu^ + r'^puft + c ^ r^n^^ 
+ j^r^{p^)turtt + jmr^-\p^)tutt + c j^r^fl 
where we have used Lemma |2. II and the following inequality found in [3]: 

{r'^ul + mr"'-'^u^){r,t)dr 



sup |«(r, t)\ < c 

a<r<b 



and the one-dimensional Sobolev inequality: 



sup \ip{r)\ < c\\ip\\h^(^i}- 

a<r<b 



forn(a,t) = 0, (2.3) 



(2.4) 
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Thus 



< 



- j^r'^ptututt - j^r'^iptuur - ptf)uu + j^r"'{p^)tUrtt + j^mr"'-^{p''')tUtt 



d 

dt Ji 
1 



r™ptn? + r^{ptuur - ptf)ut - r^{p^)tUrt - mT^-\p%ut 



+ ^ y r"'pttUt+ J r^'ipttUUr + PtUtUr + PtUUrt - Pttf - Ptft)ut 

r'^ip^'hurt- I mT'^~\p'<)ttUt + c I r^'ul^ + c [ r^f? + c 



d 

< 

dt 

1 



-r"'ptui + r^'iptuur - ptf)ut - r^{p^)tUrt - mr"'-\p'r)tUt 



- {r^pu)rtut + c\\utU^ / r"^{pf^+u'ut)+c\\ut\\U / r'^ipt+K) 



+c\\uth 



+c j^r^'lip^U^ + cj{r"'ult + mr"''\l) +cj^r"^f? + c, 
where we have used (|1.4p i . (j2.3p . Cauchy inequahty, Holder inequahty and Lemma |2.1[ This 
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together with (j2.3p . (j2.4p . Cauchy inequahty, integration by parts and Lemma |2 . 1 1 imphes 
1 



(r + mr ) 



2 7/ ^"«+2diA 



< 



1 



+ / r"^{put + ptu)utUrt + C 



d 

< , 

dt J I 



-r^'^ptui + r^{ptuur - ptf)ut - r^{p^)tUrt - mr^-\p'')tUt 



+c J r'^ip'ut + pfu'ui) + c J {r'^uit + mr^'-'uf) 
+cj^r^{\{p^)u\''+pl + fl) + c 



d 

< 

dt .ij 



^r^ptui + T^{ptuur - ptf)ut - r^{p^)turt - mr^-\p')tut 



+c\\nt\\U ( l^r^pnU J^r"^p't ) + c / {r^ut^ + mr"^--'ut) 



+cj^r^{\{p'r)u\' + pl + f^) + c 



< 



dt J J 



r'^^ptul + r^{ptuur - ptf)ut - r^{p'')tUrt - mT^'^-\p^)tUt 



I r'^ipl + f?) + l] +cj^r^\\{p'')u? + pI + fl) + c. 



Integrating the above inequality over (0, t), we have 



< 



1 



^r^ptut + r'^iptuur - ptf)ut - r^{p^)tUrt - mr^-\p^)tUt 



1 



l^r-^{pi+m+i 



+ c 



+c /* / {r^ul, + mr^-\l) [ r^{pl + ff) + c, 
JO Ji Ji 

where we have used ()1.4p i . (jl.Sp . (j2.3p . Lemma 12.11 Holder inequality and the following 

equalities: 

/ (r"^n2, + mr— 2^2)1^=0 = c [ |Vut 
Ji Jn 



itP|t=o, 



Ut\t=o = po ^ {Luq + pof (0) - VP(po) - PqUq ■ Vuo) 
= gi - uo • Vuo G -Do- 



We apply integration by parts, (j2.3p . ()2.4p . Holder inequality, Cauchy inequality and Lemma 
l2.1l to obtain 




^/ 



< 



r pUUtUrt + c 
t 



,m-2„,2\ 



J/ J/ 

1 , ,1 

m„.2 , _ m-2„.2N 



< c ( / r™/>n? 



< 



Therefore, 



JO J/ J/ 



10 J I ^ J I 

JO J/ J/ 




By Gronwall inequality, we get 



r pUtt + [r u^t + r Ut) < c. 
Qt JI 



The proof of Lemma 12.21 is completed. 
Lemma 2.3 For any < t < T , it holds 



□ 



Proof. Differentiating ()2.ip with respect to r, we get 

Urrr = Pr^t + pUrt + PrUUr + puf. + pUU^r + {p'^)rr — 



mUrr , 2m{rur -u) . /o c;^ 
^ o Pr / - P/r- (2.5) 



By ([23]), (fOj) . ([23]) and Lemma EH we obtain 

< cjir'^ul^ + mr^-Wt) + c j ^r'^\{p^)rr? +cj^r^{f + fi) + c 



^m^2 



< C. 



The proof of Lemma 12.31 is completed. 
Lemma 2.4 For any < t <T, it holds 

r"'{pl + \ip')u\') + 



□ 



r'^ult < c. 



Qt 
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Proof. By ([221), (1131), ([23]), Lemma O and Lemma [121 we get 

/" r™n2,i ^ c /" r"p2^2 ^ ^ f r^^p^ul + c [ r'^pfu'^u^, +c [ r'^p'^ulul 

Jqt Jqt Jqt Jqt Jqt 



Qt JQt JQt 



+c[ r-{plf + p'fl) 
JQt 

< C. 

Multiplying (|1.4p i by jp'^~^,we get 

{p^)t + Ip'^Ur + {p'^)rU + mjr~^p'^u = 0. (2.6) 
priDi. (I13D, (EaD, (USD, Lemma O and Lemma imply 

r'^{pl + \{p'h\')<c. 



This proves Lemma |2.4[ □ 
To sum up, we get 

Up,p')\\h? + \\{pu{p')t)\\m^ + IIU,(/'^)*t)llLp <c, P>^-, (2.7) 

and 

/ r'"(u2^ + r-'^uf + ul„ + ul^ + ti^ + r'^n^) + / r'"(pn| + n^^J < c, (2.8) 
J I JQt 

By (I12D and (USD, we get 

l|p||//2{7) + \\{Pt,Ut)\\HHl) + \\ptt\\L-^[I) + \\u\\H-i(I) + I {\\ut\\H^[I) + \\uu\\l'^{I)) < c{b,6), 

p > -• 

c 

This proves Theorem l2.1[ □ 
Proof of Theorem ll.lt 

Let b < oo, and denote Pq = Po + for 5 G (0, 1), we have 

pi^po, in H\I), (2.9) 

(/,^)^ ^ pI in (2.10) 
Let Uq be the unique solution to the equation: 



^or + ^ ) - iUn- = Pi[9i - /(O)], in /, (2.11) 



for Mob/ = 0- P-SP implies 



(no. + ^). - [(/>o)1. = Po[9i - /(O)], in /, (2.12) 
r 
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for uoldi = 0. 

From (j2.9p - (|2.12p and the standard elliptic estimates, we obtain 

uf^ uo, in H^il), as 5^0. (2.13) 
Consider (jl.4p - ()1.6p with initial-boundary data replaced by 

ip',u%=o = ipi4), in /, 

and 

u^ldl = 0, for t > 0. 

Then we get a unique solution {p^,u^) for each (5 > by Theorem 12.11 
It follows from (EZl) and (ESI) that 



\\ip\ip'r)\\HHi) + \\ipt{ip'r)t)\\mii) + \\{pL{ip'rh)\\LHi) < c, (2.14) 

and 

J ^ I ll„.<5|| I f fJl„.S |2 , L.5 



/ > -, HWhHD + Wu'WhHD + / (/Kr + Krtn < c{b). (2.15) 

Based on the estimates in (|2.14p and (|2.15p . we get a solution {p, u) to (jl.4p - (|1.6p after taking 
the limit (5 — )• (take the subsequence if necessary), satisfying 

(p,p^) G L°°{[0,T];H\I)), {pt,{P^)t) G L°°{[0,T]; H\I)), 

{pu,{p^)tt)eL--{[0,T];L\l)), p>0, ^uu e L^{[0,T]; L^{I)), (2.16) 
uG L~([0,r];F3(/)ni?i(/)), Ut€L°-{[0,T];H^{I))nLH[0,T];H^{I)). 

Since u G L°°{[0,T]; H^{I)) and ut G L°°([0, T]; F(J(/)), then we get u G C([0, T]; //^(j)) 
(refer to ffj). By (ll.4j) i . (12. 6j) and similar arguments as that in [HEj, we get 

pC^C{[0,T];HHl)), pteC{%T]-H\l)), (2.17) 

and 

p^ G C7([0,T];F2(/)), {p^), G C([0,T];Fi(I)). (2.18) 
Denote G = {ur + ^ - p'<)r + pf. By dH]) and (IZTHD . we have 

G = put + puur G L2([0,r];F2(/)), 

Gt = {put + p^x^x,)t G l2([0, T];L\I)). 

By the embedding theorem (ilj), we have G G C([0, T]; i7i(/)). Since pf G C([0, T]; i7i(/)), 
we get 

(ur + ^-p'')^GC{[0,T];H\l)). 



This means 



+ — -p^ G C7([0,r];i72(/)). (2.19) 
r 
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By (I238D and (I2J9D . we get 

Ur+'^e C{[0,T];H''). 
This together with u e C([0, T]; impHes 

u(^C{[Q,T\]H^{I)). (2.20) 

([LiD o. (l2Tril . (f2J8]l and (l2:20D give 

(HteC([0,r];/?^(/)). (2.21) 

Denote 

p(x, t) = p{r, t), u(x, t) = u{r, t)-. (2.22) 

r 

It follows from (l2TB - (l2J8l ) and (l2:20]l - (12:22]) . we complete the proof of Theorem O for 
6 < oo. For 6 = oo, we can use the similar methods as that in [3] together with the estimates 
()2.7p and ()2.8p uniform for 5 to get it. We omit details here for simplicity. □ 

3 Proof of Theorem D 

Similarly to the proof of Theorem 11.11 we need the following auxiliary theorem. 

Theorem 3.1 Consider the same assumptions as in Theorem M.SX and in addition assume 
that Pq > 5 > Q and b < oo. Then for any T > 0, there exists a unique global classical 
solution {p,u) to |i.^[ j-( [T7^) satisfying 

p€C{[0,T];H^iI)), p>i uGCi[0,T]-H'^iI)nH^iI))nL\[0,T];H^I)), 
uteC{[0,T];H^I)n H] (!) ) n L^{[0,T]- H^I)), 
utt € C{[0,nHUl))nL^{[0,T];H^I)), uut G L^{[0,T]; L\I)). 

Similarly to the proof of Theorem 12. H Theorem 13.11 can be proved by some a priori estimates 
globally in time. Since ()2.7p and ()2.8p are also valid here, we need some other a priori 
estimates about higher order derivatives of {p, u). The generic positive constant c may depend 
on the initial data presented in Theorem 1 1 . 2 1 and other known constants but independent of 
6 and b. 

Lemma 3.1 For any < t < T , it holds 

I r^[pl„. + pI, + \{p^\rr? + W)rrt?] + [ r^[plu + \{p')rU? + ul„] < C. 
J I JQt 

Proof. Taking derivative of order three on both sides of (jl.4p i with respect to r, we have 

_i o -1 o -1 3mpUrr QmpUr _i 
p^^^t = —mr pUj-rr — 'imr PrUrr — oTTir PrrUr H k n rnr PrrrU 

3mprrU QmprUr GmprU 6mpu 
H 7^ ~h n n H" "a Pfj'rpj'U ApfffZif ^Pfj-iiff 
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Multiplying (j3.ip by r"^prrr, integrating the resulting equation over /, and using integration 
by parts, we have 



6mpur _i Smprr-u GmprUr GmprU 6mpu 
o mr PrrrU H 5 1 5 ^ 1 j 4/3 



—AprU 



+ \ j r'^Prrr'^r + \ j "^'^"^ Vrrr'" " j r"^ PPrrrUrrrr ■ 



By Sobolev inequality, (j2.3p . ()2.4p (j2.7p . (|2.8p and Cauchy inequality, we get 

r'^Prrr < ^ Prrr + C ^ r'^U^rrrr + C (3.2) 



Similarly to (j3.2p . we get from (j2.6 



^r™|(/5^)„,|^ <c j^r^\{p^)rrrV + c j^r^ui,,, + c. (3.3) 



(it 



By ([321) and ([331), we have 

r™(p2^, + |(/9^)„.P) < cj^r^'ipl, + \{p^)„r?)+c j^r^ul,,, + c. (3.4) 
Differentiating (12. 5p with respect to r, we have 

= PrrUt + 2prUrt + /OUrrt + {prUUr + /OU^ + pUUrr)r + (p''^) 

+3mr~'^Urr — Gmr'^Ur H ^ Prr/ — 2/9r/r " P/rr. (3.5) 

From (133]), dOD, dMD, and ([231), we obtain 

j^r'^ul,„<c j^r'^\{p'')rrr? + c j^r'^pul.t + c j^r'^fl (3.6) 
By dH, dMl), (Eai) and dlZ]), we get 

■J^r^iplr + l(/'^)rrrP) < cj^T^{pl„ + |(/>^)„.|2) + C^r^tX^,, + C^r"7'r + C 

By Gronwall inequality and (|2.8p . we obtain 

^r^{pl„ + \{p')rrr?)<C. (3.7) 

It follows from ([Ll])i, ([23]), (HID, ([SSD, ([22]), ([M]), dSZI) and Lemma O that 

/ r'^iplrt + l(/^")r.iP] + f r^[plu + \{P''U? + n2_] < c. 

The proof of Lemma |3. H is completed. □ 
Lemma 3.2 For any T > 0, we have 

\\{\/^)r\\L^{QT) + \\{y/p)t\\L^(QT) < C. 
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dt 



Proof. Multiplying (ll.4p i by , we have 



TTIT ^ 1 
(\/p)t + -^VPU + -^/pUr + {./p)rU = 0. 



Differentiating (|3.8p with respect to r, we get 



iVPjrt + -^{VP)rU + -^r—^/pUr + -{^/p)rUr + -^fP^rr + \^JP) 



mr 



ruy/pu 3 



2r2 ' 2^^"^'"' ' 2 



Denote h = {^)r-, we have 



ht + hrU + h[ —U + -M^) H —^/pUr h -^/pUrr = 0, 



TTLy/pU 1 



which implies 



/lexp 



+ 



X exp 



2 2r2 ' 2 



/o 



2 + ^ ) (?^(T,y),T)(ir 



where r{t,y) satisfies 



^=u{rit,y),t), 0<t<s, 
jis,y) = y. 
Integrating (|3.9p over (0,s), we get 



h{y, s) = exp 



+ (r(r,y),r)dr ) Mr(0,y),0) 



2 2r2 ' 2 

* / mr~^u 3u,. , , , , , , 
xexp( / ( — h — ) (r(T,y),r)dr 



V 2 2 
This together with ([23]), ([23]), ([IZD and imphes 

ll(\/p)r||L->(Q^) < C. 

From dSSI), (I3J0D . (lOD . (ITiD . (|2TD and (gSD, we get 

ll(\/p)t||L-'(QT) < C. 

The proof of Lemma 13.21 is completed. 



dt. 



Lemma 3.3 For any < t < T, it holds 

/ + P'^rrt + Urrrr) + / ^^{P^'^ttt + '^rrrt) — ^■ 

J I JQt 
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Proof. Differentiating ()2.2|) witli respect to t, we get 

PUttt + PttUt + ^ptUtt + PttUUr + ^ptUtUr + 2ptUUrt + pUuUr + 2pUtUrt 

+pUUrtt + (/5'^)r« = ^^rrtt + mr~'^{rUrtt " + + '^Ptft + /^/tt- (3.11) 

Multiplying (I3.1ip by r^p^utt, integrating the resulting equation over /, and using integration 
by parts and Cauchy inequality, we have 

1 d 



2dt 
1 

~2 



m 3 2 , / m 2 2 , m~2 2 2 \ 

r p Utt+ [r p u^tt + mr p Utt) 
I Jl 

T^p'ptul - [ r^p^u 
I J I 



PttUt + PttUUr + 2ptUtUr + 2ptUUrt + pUttUr 
+2pUtUrt + pUUrtt + {P"')rtt -2 j r"" pPrUftUrtt + j r"^ p'^Uttifttf + 2ptft + pftt) 

< cj^r^pul + ^l^r"^p\%, +cj^r^WU\^ - ^ j^r^ pUrttVpuu{^)r 
+cf^r"^{pltP + p'tf^ + p'fl) + c 

< c r^pul + r^p'ul,, + c r^\{p<)rtt? + r^^p'u^tt + ^ r"^fu + c, 

where we have used (|2.3p , ()2.4p , ()2.7p , ()2.8p , Lemma 13.21 and Cauchy inequality. 
Thus, 



It 



r^p\l + / [r^p'ulu + mr^-'p^ul) 



I 



< cj^r^pul + cjn{p')rtt? + j^r^fl + c. (3.12) 



Integrating (j3.12p over (0, t), and using (|2.8p and Lemma [3JJ we get 

i-t I- I- 

,.m„3„,2/,\ I / lf^rri2^,2 . ^„m~2 2„,2\^ I ^rn J,„,2 



T^p'utM + / / {r^p'uU, + mr^-'p'ui) < / v^p'uiM + c. (3.13) 
/ JO Ji Ji 

By (gH), (gSD, (lESl) and ^V^gi G L^, we have 

I^r"^p\li0)<c. (3.14) 

(I3J3]1 and ^AM give 

/ r^p\l + / {r^p'ul, + r—V'n?, ) < c. (3.15) 
By (I22|), Cauchy inequality, (H^D, ([231), ([22]), 1^ and (f3l3D . we have 

r"'pulrt < c f r'^ppWt + c / r'^p'^ul + c [ r'^'ppiu^ul + c [ r'^p^ulul 



+cj^r"'pWlt + cj^r"'p\ip^)rt\^ + cJ^ r'^-^pult + c r^^-^pu^ 



+c J^r^pptr + c j^r^p-'ft 
< c. (3.16) 
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Differentiating (j2.2p witli respect to r, we get 

Urrrt = PrtUt + 'i'{^)r\fpUtt + PtUrt + PUrtt + Prt^^^^r + PrUtUr + Pr^ri + Pt^*r 
-\-2pUrUrt + PtUUrr + pUfUrr + pUUrrt + {p'^)rrt — mr~^Urrt + 2mr~'^Urt 



2mu 



Prtf - Ptfr - Prft - Pfrt- (3.17) 



By dnSD, (IXTTD . (lOI) . (I23D, (I22D, dlSI), Lemma O and LemmaO we have 

7* XL^^^-j^ C« 



3.6p . (j3.16p and Lemma |3. II immediately give 



r'^ul,,, < c. 



The proof of Lemma 13.31 is completed. □ 
Lemma 3.4 For any < t < T , it holds 

I ^ [Prrrr ~^ \ ip'^}rrrr\ ] ~l~ / ^ U^j.^^.j, ^ C. 
J I JQt 

Proof. Differentiating ()3.ip with respect to r, we get 

_^ _i Ampurrr \2mpu„. ISmpUr 

AmprrrU 12mprUrr 12mprrUr 2Ampi.Ur 



+4mr Prrr-Ur + mr Prrrr^i 5 ^ 5 1" 



^2 ^2 ly*^ 



12mprrU GmpUr 2AmprU 2Ampu 

H T 7 7 I c l~ pUrrrrr ~l~ oprUrrrr 

ry\j j^rt j^rt iy*KJ 

Multiplying (I3l8]) by integrating over /, and using integration by parts, (12. 3p . (12. 4p . 

([TTp . (USD, Lemma EH Lemma [331 we get 



Similarly, we have 

I ^r™|(p^)„.,,|2 < c^r™|(p^)„„|2 + cj^r^ul,,,, + c. (3.20) 
By (l3l9D and (fCTjl . we obtain 

jj^r^{pl„ + |(p^)„,„,|2) < cj^r^{pl„., + |(p^)„„|2) + cj^r^ul,,, + c. (3.21) 
Now we estimate JjT^u^^^^^. Differentiating (|3.5p with respect to r, we have 

Ut + SprrUrt + '^Pr-Urrt + pUrrrt + {prUUr + + pUUrr)rr + ip"')rrrr 

-1 < -2 -■^ < -4 2Amu 
—mr Urrrr + 4mr Uj-rr — 12mr Urr + 2Amr Ur g — 

Prrrf ^Prrfr ^Prfrr Pfrrr- (3.22) 
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It follows from ([3221), <^M, (I23I), (1221), dZH]), Lemma O and Lemma that 

I Ji Ji Ji 



+c ^ (/2 + fl + C^) + c. (3.23) 



Since 

„m„2„.2 A I m\i r-\ |2„.2 



d 
di 



I J I 

This together with ()2.7p . (j2.8p . Lemma 13.21 and Lemma 13.31 gives 

j^T^ul„, < cj^T^p\l,,^ + cj^r^\{p')rrrr? + cj^r^{f, + fl + + C (3.24) 

Substituting ([3:21) into dMO, we obtain 

r (^Pj,^j,^ -\- \ {p'^)rrrr\ ) 

< cj^r^ipl,, + |(p^)....|2) + cjT^ul„^ + cj^T^{f, + fl + /2,,) + c. 
Using Gronwall inequality and Lemma l3.3[ we get 

r^ipl,, + \{p"')rrrr\^)<C. (3.25) 

It follows from U^^ . IK25\i and Lemma [331 that 

This proves Lemma |3.4[ □ 
From ([131)1, ([231), (1231), ([221), (1221), ([231) and Lemmas [33K331 we immediately get the 
following estimate. 

Lemma 3.5 For any < t <T, it holds 

r"" [p'rU + \iP^)ru\'+plrrt + \iP')rrrt\'] 



I 



+ [ r^[plt + \{p'')tu\'+plu + KP 
JQt 

Lemma 3.6 For any < t < T , it holds 

f {r"^p\'ru + r-'-'p'ui) + [ r^p^i, < c 
Ji JQt 
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Proof. Multiplying (j3.1ip by r'^p'^uut, integrating the resulting equation over /, and using 
integration by parts, we have 



= 2 y r"^ p^ ptU^tf- - 4 j r"" p^ PrUrttUttt - j r'^p'^Uttt PttUt + 2ptUtt + PttUUr + 2ptUtUr 
+2ptUUrt + pUttUr + 2pUtUrt + pUUrtt + (p^)rit j '^'^"'''^ Pf^tt 

+ 1^ r^p\tu (pttf + 2ptft + pftt) 

< c r^p'ul, - 8 r^piuutp{^p)rUru + + ^ + c ^ r™(/2 + + /2 ) 

< c T^p'ul, + r^p^ul, + cj^r^{f + f? + ) + c, 

where we have used ()2.3p . (|2.4p . (|2.7p . (|2.8p . Lemma [3121 Lemma [331 Lemma [331 and Cauchy 
inequality. 
Thus, 

I r'^p'nl, + I y (r-^/n?,, + mr- < c ^ r^p^^^, + cjj^{f + f? + /.^J + c. 
By dLl]), ([221) and poV^gi G we have 

m 5 2 I /■ m 4 2 , m~2 A 2 \ ^ 

r P + [r P Urtt + r p Uu) < c. 
J I 

The proof of Lemma 13.61 is completed. □ 
Lemma 3.7 For any < t <T, it holds 

J/ JQt 
Proof. From (fXTTp . (fOD . (HZ]), (US]), Lemmas EIES] and Lemma ESI we get 

< c. 

This combining (j3.24p and Lemma 13.41 gives 

r"^|a^up < c. 

By (pmi) . (1^ . (ESD, (gTI, (I2S1), Lemma EH Lemma E3] and Lemma EE we get 



/ r^P^rtt < c[ r-(/2 + /2)+c 



IQt JQt 

< c. (3.26) 
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Differentiating (j3.17p witli respect to r, we liave 

Urrrrt = PrrtUt + '^PrtUrt + PrrUtt + ^PrUrtt + PtUrrt + PUrrtt + PrfUUr + p^UtUr 
+PrUUrf + Pfuf. + 2pUrUrt + PtUU„ + pUtU„ + pUU„-t + {p'^)rrt 

-mr~^Urrt + 2mr~'^Urt - ^"Y** - Prrtf - '^Prtfr " Prrft " Ptfrr (3.27) 

J r 

— 2prfrt — Pfrrt- 

By (lOHD . (lOTD . (ITaD . (ITID . (12771) . dZBD, Lemmas EJESI and LemmaESl we obtain 

r'^pulrrt < C r"'p\^)r\^ultt + C I r'^P^U^.rtt 

Qt JQt JQt 



-c [ r^{fl + fl + f„^)+c 

JQt 



< C. 

The proof of Lemma 13.71 is completed. □ 
Lemma 3.8 For any < t < T , it holds 

r™(|a,5/9p + \d^.{p^)\^) + [ r'^ldSl^ < c. 
I JQt 

Proof. Differentiating ()3.18p with respect to r, we obtain 

d^Pt + mr"^ pUrrrrr + rnr~^ p^Urrrr ~ mr""^ pu„rr + (^4:mr~^ PrUrrr + ()mr~^ PrrUrr 

AmpUrrr 12mpUrr ISmpUr , _i \ _i _i 
^ 1 5 -. h 4mr PrrrUr + mr PrrrrrU + mr PrrrrUr 

ryZ ipO rp^ J ^ 

5mprrrrU AmprrrUr SmprrrU 
^2 



12mprUrr 12mprrUr 2AmprU' 



+ 



j^2 



12mprrU 6mpUr 24mprU 2Ampu 



+ 



j^3 IJ-*^ IJ-*^ rjr*^ 



Urr + ea^pUr + ud'^P = 0. (3.28) 



Multiplying (fM]) by r™9^/), and using ([23]), ([231), (1221), (ESI), Lemma EH Lemma [33 
Lemma 13.41 and Lemma 13.71 we have 



J^J^r^'ldrpf < c j^r"'\dlp\^ + c j^r'^\dlu\^ + j^r'^{\dlpWu + c 
< cj^r"-'\d^p\'^ + cj^r'^\d^u\'^ + c. 



(3.29) 



Similarly, we get 

d 
di 



j^r'^\dl{p^)\'^ <c j^r'^\dl{p^)\'^ + c j^r'^\d%\'^ + c. (3.30) 
Differentiating (|3.22|) with respect to r, we get 

til — PtTVT I ^Prpiprp Uji^-f^ I Q Pfi^ 'Ufi^f^ j 4jf?^ '^TTT't I pXiif'ff'ff^ I (^Pi^ tiiUif I pUjrp I P'WUfY' ^ TVV 

—mr~^Urrrr + 4mr~'^ii^^^ — 12mr~'^n^^ + 24mr~^u^ = — I 

Prrrrf '^Prrrfr ^Prrfrr ^Prfrrr Pfrrrr- (3.31) 
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(lOTD . don . (HaD, (I22D, d^H]), Lemma EH Lemma ESI Lemma EH and Lemma E21 imply 
J I J I J I J I J I 

+C j r {fj.j. + fj-rr ~^ frrrr) ~^ ^- (3.32) 

It follows from (fCTjl . (^150]) . (pTTSD . (ITSD . Lemma E31 Lemma EZ3 and Gronwall inequality 
that 

J^r^[\d',p\' + \dUp^)\']<c. (3.33) 
From (|3.32|) . (|3.33|) . (|2.8|) . Lemma E31 and Lemma E21 we obtain 

[ r'^ld^ul^ < c. 
Jqt 

The proof of Lemma 13.81 is completed. □ 
It follows from l^i, (EJD, (EZl), (USD and Lemmas EID that 

\\ip,p'')\\H^ + WiPU ip^)t)\\Hf + mVP)r, iVP)t)\\L-iQT) <C, P>^-, (3.34) 

and 

/mi- 3 2 I 4 2 , -2 4 2, 2 2 , 2 , 2 , — 2 2 , i o5 |2 , i o4 |2 
r [p + p u^tt + r p + p u„^t + pu^^t + u^t + r + \d^u\ + \d^u\ 

Wrrr + ^Ir + ^^r + r^\^) + / r™(pUj2t + p^ul^^ + /J^Uj^ti + p^U^rtt + P^lrTTt 

Jqt 

+U^rt + ^rrrt + l^r^l^) ^ C. (3.35) 

By (lOil) and dOS]) . we get 

\\p\\h'o{i) + \\pt\\m{i) + \\u\\hhi) + \\ut\\m{i) + \\utt\\m{i) + {\\u\\h6(i) + htWH^i) 
+ \\uu\\h2^j) + \\utu\\l2^j)) < c{b,5). 
This proves Theorem 13.11 □ 

Proof of Theorem 11.21 

Since ()3.34p and ()3.35p are uniform for b and 5, it suffices to prove Theorem 11.21 for the 
case b < oo. We follow the strategy as the proof of Theorem 1 1 . 1 1 and use Theorem 13. 11 After 
taking 5 — )■ (take subsequence if necessary), we get a solution {p,u) of (11.4|) - (11.6p satisfying 

(p,p^) G L-([0,T];F5(/)), {{^)r,{VP)t) e L^{Qt), 
{pu{p<)t)eL^{%T]-H\l)), u€L--{[0,T];H'{I))nLH[0,T];H^I)), 
ute L^{%T]-Hl{I))rM\%T]-H'\I)), 

{^d^ut,pd?ut) e L^{[0,T];L^{I))nL\[0,T];H\l)), (3.36) 
^d^ut G L^{Qt), P^utt G L^{Qt), P^uttt G L^iQr), 
pluu€L^i[0,T];L^{I))nL^{[0,T];HUl)nHHl)), 
p'uueL--{[0,T];H^,{I))nL\[0,T];HHl)). 
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It follows from u G L'^{[0,T]; H^{I)) n L'^{[0,T]; H^{I)), ut G L\[0,T]; H'^{I)), M 
and the similar arguments as [Sj E] that 

Denote p(x, t) = p{r,t), u(x,t) = u{r,t)^, then (p, u) is the unique solution to (ll.ip - (jl.3p 
with the regularities in Theorem 11.21 The proof of Theorem 11.21 is completed. □ 

Remark 3.1 By our method, it seems that the regularities of u could not be improved to 
L°°{[0,T];D^{n)) and L'^{[0,T]; D^{n)), even if the systems ([ri])-([L3D have initial data, f 
and gi smooth enough. More precisely, based on ()3.34p and ()3.35p . the next two a priori 
estimates about u for (jl.4p - p.6p are 

/ r'^p'ult + [ r^p\ul,u + r-\L) < c, (3.37) 
J I JQt 

and 

/ r-p8(^2^^^ + ,-2^2^) + f r^p^ulu < c. (3.38) 
Jl JQt 

p.37p and (|3.38p respectively implies 

/ r'^p\dtu\^ + / r^iaj-up < c, (3.39) 
J I JQt 

and 

[r'^p^\dlu\^+ [ r'^p\d^uf<c. (3.40) 
Jl JQt 

Just because of the appearance of the vacuum, we can't obtain the regularities u in L°°(0, T; D^{^)) 
and L'^{id,T;D^{Vl)). 
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